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ABSTRACT 


The  power  density  spectrum  of  the  received  acoustic  waveform  after  reflection  off  the  time-varying 
random  sea  surface  is  evaluated  at  an  arbitrary  observation  point  in  the  farfield.  For  a  monochromatic 
transmitted  signal  and  a  narrow-band  Gaussian  surface -height  variation,  the  received  acoustic  spectrum 
is  shown  to  consist  of  an  impulse  at  the  transmitted  acoustic  frequency  plus  sidelobes  centered  at  fre¬ 
quencies  separated  from  the  transmitted  fvquency  by  multiples  of  the  surface  center  frequency.  The 
powers  in  the  coherent  component  and  scattered  sideband  components  of  the  received  pressure  waveform 
are  evaluated  in  terms  of  the  surface  roughness  and  spatial-temporal  correlation  function  of  the  surface. 
For  the  special  case  of  elliptical  contours  of  iso-correlationat  zero  time  delay,  the  sideband  powers  and 
scattering  strengths  are  evaluated  in  terms  of  two  fundamental  parameters  that  include  the  geometry  of 
the  experiment,  the  incident  acoustic  frequency,  the  root  mean  square  (rms)  surface  height,  and  the 
surface  correlation  distances. 

The  rms  bandwidths  of  the  sideband  scatter  components  are  evaluated  for  small  surface  roughness  and 
shown  to  be  approximately  proportional  to  the  square  root  of  the  sideband  number.  Numerous  examples  of 
sideband  scattering  strengths  for  a  variety  of  spatial  correlation  functions,  including  exponential  and 
Gaussian  decay  as  special  cases,  are  given. 
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SPECTRUM  OF  A  SIGNAL  REFLECTED  FROM  A 
TIME -VARYING  RANDOM  SURFACE 


1.  INTRODUCTION 

The  purpose  of  this  study  is  to  determine  the  spectrum  of  an  acoustic 
signal  reflected  from  a  time-varying  random  surface.  Of  the  many 
scientists  who  have  investigated  scattering  from  such  a  surface,  we 
will  mention  only  a  few.  Rayleigh  considered  a  sinusoidal  fixed  surface 
andsolvedfor  the  reflectedpressure  for  a  given  incident  monochromatic 
wave.  Eckart's  classic  work  £lj  on  scattering  from  a  random  surface 
employed  the  Kirchoff  method;  i.  e. ,  he  assumed  that  the  surface  could 
be  represented  by  the  tangent  plane  to  the  surface  at  each  point 
of  incidence.  Starting  with  the  Helmholtz  equation,  Beckmann  and 
Spizzichino  [2]  also  used  the  Kirchoff  approximation,  but,  in  their  deri¬ 
vation,  they  differentiated  with  respect  to  the  normal  to  ihe  surface 
instead  of  the  vertical  axis,  andused  an  integration  by  parts  to  eliminate 
edge  effects.  Thus,  they  obtained  a  different  geometric  factor  from  that 
of  Eckart  [l]  for  the  reflected  pressure.  LaCasce  and  Tamarkin  [3]  per¬ 
formed  experimental  studies  on  the  reflection  off  a  fixed  sinusoidal 
pressure-release  surface.  Horton  and  Muir  [4]  also  performed  experi¬ 
mental  studies  on  the  scattering  from  a  rough  random  surface,  compared 
their  results  with  a  modified  version  of  the  Eckart  theory,  and  found 
good  agreement.  Wagner  [5]  considered  a  Gaussian  surface,  andusing 
the  parameters  of  root-mean-square  (rms)  height  and  correlation  of 
surface  heights,  and  treating  the  incident  waves  as  rays,  determined 
the  amount  of  shadowing  as  a  function  of  the  grazing  angle  of  the  incident 
acoustic  wave.  From  these  results,  a  correction  value  for  the  shadowing 
effect  was  made  by  Wagner  [5]  to  the  Eckart  [l]  and  Beckmann  and 
Spizzichino  [2^  theories.  Marsh  [6]  also  utilized  the  Rayleigh  method 
and  considered  the  signal  reflected  from  a  random  surface  to  consist 
of  an  infinite  set  of  plane  waves.  By  applying  boundary  conditions  on 
the  surface,  he  solved  for  the  scattering  coefficient  of  the  surface 
reflection. 

The  number  of  investigations  of  a  time-varying  random  surface  is 
considerably  less  than  that  of  a  fixed  surface.  Of  these,  we  will  mention 
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only  the  works  of  Roderick  and  Cron  [7]  and  Parkins  [8]  .  Roderick 
and  Cron  considered  a  sinusoidal  traveling  surface  wave,  where,  for  a 
single-frequency  incident  acoustic  wave,  the  frequency  of  the  reflection 
in  the  specular  direction  is  the  same  as  the  incident  frequency,  and  the 
frequencies  in  the  nonspecular  directions  are  equal  to  the  incident  fre¬ 
quency  plus  multiples  of  the  surface  frequency.  They  compared  the 
theoretical  predictions  of  the  amplitudes  of  these  components  with  ex¬ 
periment.  Parkins  considered  a  time-varying  random  surface,  a  joint 
Gaussian  distribution  of  heights,  anda  Neumann-Pierson  surface -height 
spectrum.  For  both  small  and  large  surface  roughness,  he  obtained 
the  scattered  acoustic  spectrum. 


In  our  work,  we  start  with  the  assumptions  used  b> 


•  theKirchoff  method  [2,  p.  20]  _  the  methodof  physical  optics  — 
is  used, 

•  the  insonified  area  is  in  the  farfield  (Fraunhofer  region)  of  the 
source, 

•  the  observation  point  is  in  the  farfield  of  the  insonified  area, 

•  the  surface  is  pressure-release, 

•  no  shadowing  or  multiple  reflections  occur, 

•  a  directional  source  is  used,  and 

•  the  surface  particle  velocity  is  small  compared  with  the  speed 
of  sound. 

Starting  with  the  Helmholtz  integral  and  a  given  incident  single¬ 
frequency  acoustic  source,  we  obtain  the  (complex)  pressure  at  the 
observation  point.  The  surface  height  is  designated  a  function  of  position 
x,  y,  and  time  t;  the  first-  and  second-order  probability  distributions 
of  height  are  assumed  independent  of  absolute  position  and  time  and 
dependent  only  on  tie  differences  in  the  x  and  y  coordinates  and  time. 
That  is,  the  surface  is  assumed  homogeneous  and  stationary.  The 
autocorrelation  function  of  the  received  acoustic  complex  pressure 
is  obtained  by  using  these  properties.  The  Fourier  transform  of  the 
autocorrelation  function  then  yields  the  received  acoustic  power  spectrum. 


The  equation  for  the  power  spectrum  is  then  specialized  to  the  case 
of  a  joint  Gaussian  distribution  of  surface  heights  and  a  narrow -band, 
surface-height  spectrumat  a  point;  for  this  case,  the  spatial-temporal 
correlation  function  of  the  surface  can  be  represented  as  a  sinusoidal 
oscillation  in  time-delay  with  a  slowly  varying  amplitude  and  phase.  It 
is  then  assumed  that  the  iso-values  of  the  surface -correlation  function 
is  an  ellipse;  for  this  case,  it  is  shown  that  the  received  acoustic  spec¬ 
trum  consists  of  a  series  of  spectral  lobes  or  sidebands,  each  separated 
from  the  incident  acoustic  frequency  by  multiples  of  the  center  fre¬ 
quency  of  the  surface  variation,  plus  an  impulse  at  the  incident  acoustic 
frequency.  The  relative  powers  in  these  sideband  components  are 
evaluated  in  terms  of  two  fundamental  parameters:  the  first  is  the 
surface-roughness  parameter,  and  the  second  is  related  to  the  corre¬ 
lation  distances  of  the  surface  (the  horizontal  separation  for  which  the 
surface  correlation  falls  to  e-1  of  its  origin  value).  The  powers  in 
the  sidebands  are  numerically  evaluated  and  plotted  versus  these  two 
parameters,  for  both  specular  and nonspecular  directions,  for  a  variety 
of  spatial  correlation  functions,  including  exponential  and  Gaussian 
decay  as  special  cases. 


2.  REFLECTION  FROM  A  GENERAL  TIME-VARYING  SURFACE 

2.  1  RECEIVED  PRESSURE  WAVEFORM  REFLECTED  FROM 
A  TIME-VARYING  SURFACE 

The  geometry  of  the  scattering  experiment  is  depicted  in  Fig.  1. 


Pig.  1  -  Scattering  Geometry 
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The  source  o  1  acoustic  radiation  is  at  Q,  and  the  observation  point  is 
at  A.  The  svmbol  O  is  an  origin  of  coordinates.  Let  the  unit  vectors 

4  S  —A’  ^ 

in  the  x,  y,  z  directions  be  i,j,k,  respectively,  and  let  uQ  and  uA 
be  unit  vectors  from  O  in  the  directions  of  source  Q  and  receiver 
A,  respectively.  Then, 

“Q=V  +  bQ?  +  cQt  » 

-»  T*  . 

UA  =aA*  ^bA>  +CAk  ’ 


where  (a  ,  b  ,  c  )  and  (a  ,  b  ,  c  )  are  the  direction  cosines  of  the 

0  Q  Q  AAA 

source  and  receiver,  respectively. 

Let  this  sums  of  direction  cosines  be  denoted  by 


a  =  aQ  +  *A 
b  *  bQ  +  bA 


C=CQ  +  CA 


(2) 


The  specular  direction  corresponds  to  a  -  b  =  0;  i.  e. ,  aA=  -aQ, 
bA  =  -bQ,  and  cA  =  cQ.  It  is  the  direction  in  which  all  the  reflected 
energy  would  occur  for  a  mirror-like  surface. 


The  height  of  the  reflecting  surface  at  position  x,  y  fluctuates  with 
time  and  is  denoted  by  £ (x,  y,  t).  The  average  surface  height  corresponds 
to  £  =  0.  Under  the  assumptions  given  in  Section  1,  the  (complex)  re- 
ceivedpressure  at  A,  for  a  single-frequency  excitation  exp(i2irf  ,t)  at 
Q,  is  given  by  [8,  Eqs.  (1)  and  ( 10) 1 ;  1,  Etjs.  (1)  and  (6)2;  4, 

9,  Eqs.  (6.9)  and  (6.24)] 

B 

p(t)  =  -i  e*p[i(2pf(lt  -k8Ro)l - ffdx dy  Pj  (s,y)  • 

\,Ro 

e*ptik,  (ax  +  by  +  c  £(x,  y,  t)  )  ]  , 


Eq.  (9); 

(3) 


where3 

f  =  acoustic  frequency, 

=  acoustic  wavelength, 

lThe  factor  <f>+  is  missing  in  (10)  of  Ref.  8. 

2The  exponent  inside  the  integral  in  (6)  of  Ref.  1  should  be  negative;  see  also  Ref.  4,  footnote  6. 

^Strictly,  the  argument  t  in  £  should  be  replaced  by  the  retarded  time  taken  by  a  signal  to 
travel  from  a  reflecting  point  (x,y)  to  A.  However,  for  a  slowly  fluctuating  surface  that  does  not 
change  much  in  the  litre  taker,  for  sound  to  ptopagate  across  the  illuminated  region,  a  good  approxi¬ 
mation  is  t-Rc/v,  where  v  is  the  propagation  velocity.  The  delay  R0/v  has  been  dropped  in  £  in 
(?)  for  notational  simplicity.  Also  integrals  without  limits  are  over  the  range  of  nonzero  integrand. 


4 


ka  =  acoustic  wave  number  (2 ff/A  ), 

Ro  -  distance  from  origin  O  to  observation  point  A, 

^i(x.y)  -  spatially  dependent  component  of  the  incident  pressure  field  on  the  reflecting 
surface  (without  the  phase  factor  due  to  propagation  over  the  distance  from  Q 
to  O),  and 

B  -  real  scale  factor  depending  only  on  the  geometry  of  the  experiment. (See  Ap¬ 
pendix  A  for  further  details  on  this  factor.) 

(Since  some  authors  assume  an  excitation  of  the  form  exp(i2trf  at), 
differences  in  sign  with  (3)  will  occur  in  those  references). 


For  a  general  signal  s(t)  transmitted,  with  voltage  density  spectrum 
S(f),  the  received  pressure  waveform  r(t)  is  given  by 


r(t)  =  /df a  p(t;f4)  S(£a) 

=  /dfaexp(i2nf#t)  H(f#;t)  S (f a ) 


(4a) 


where 


B 

H(fa;t)  =  -i  exp(-ikaRQ)  - ffdxdy  p.(x  y)  • 

\Ro 

exp[ilct(ax  +  by  +  cC(x,y,t))]  . 


(4b) 


We  have  indicated  explicitly  the  dependence  of  receivedpressure  p 
on  acoustic  frequency  fa(p.  is  also  dependent  on  f a )  and  used  (3).  The 
representation  in  (4a)  shows  that  the  surface  can  be  viewed  as  a  linear 
time-varying  filter  on  transmitted  acoustic  waveforms  with  instantaneous 
transfer  function  H(fa ;  t)  given  by  (4b). 


For  a  real  sinusoidal  acoustic  signal  transmitted, 

s(t)  -  cos(2n-fat)  =  Re  {  exp(  i2rrfat)l  ,  (5) 

the  real  receivedpressure  is  given  by  Relp(t)!  .  The  power  in  the  real 

received  signals  is  ______  __________  _ 

Re2(p(t)(  =  i[p(t)  +  p* (t)  1 2  =~jp(t)|  2  ,  ^ 

where  we  have  used  the  fact  that  p2(t)  =  0;  this  follows  from  (3)  upon 
noting  that  if  £  is  stationary,  then  p(t)  is  stationary,  narrow-band,  and 
centered  around  frequency  fa  and  has  its  spectral  content  confined  to 
positive  frequencies.  (See  Appendix  B. )  Under  the  reasonable  assumption 
that  surface  height  £(x,  y,  t)  varies  slowly  with  time,  as  compared 
with  exp(i2irfat),  p(t)  is  narrowband.  Thus,  attention  can  be  focused  on 

4  For  example,  surface-height  variations  typically  have  spectral  content  in  the  neighborhood 
of  fractions  of  Hertz,  whereas  acoustic  frequencies  are  of  the  order  of  tens  of  Hertz  (and  greater). 
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complex  pressure  p(t)  and  appropriate  real  parts  or  factors  of  l/2 
applied  later  when  necessary. 


2.2  COHERENT  AND  SCATTER  COMPONENTS  FOR  A  RANDOM 
SURFACE 


The  received  pressure  for  a  single-frequency  excitation  is  given 
by  (3).  The  coherent  component  [2,  Section  7.3]  of  this  waveform  is 
defined  as  its  mean  value  (ensemble  average  over  all  possible  surface 
states): 

_  B 

pc(t)  S  p(t)  =  -i  exp[  i(2fffst  -k4Ro)l - jjfdx  dy  p.  ( x,y )  ■ 

\Ro  (7) 

Mp[ika(ax  +  by)]  f ^(k^c)  , 

where 

f ^ (k a c)  s  exp[ikac  £(x,y,t)] 

= /d  ^exp(ikac^)  qj(C)  (8) 


is  the  first-order  characteristic  function  (CF)  of  the  surface  wave 
height.  Since  we  are  assuming  a  homogeneous  stationary  surface,  qj 
is  the  first-order  probability  density  function  (PDF)  of  the  surface 
wave  height  and  is  independent  of  absolute  position  x,  y  and  time  t. 
The  first-order  CF  contains  all  first-order  statistical  information 
about  the  surface  wave -height  variation  since  it  is  a  Fourier  transform 
of  the  first-order  PDF. 


If  we  define  the  double  integral  in  (7)  as  P, 

//dxdy^.  (x,  y)  exp[ika(ax  +  by)]  =  P(kaa,kab)  ,  (9) 


the  coherent  component  is  given  by 

B 

pc  (t)  =  -i  exp[i(2fff#t -kaRo)] - P(kaa, kab)  f ^  (kac)  *  (10) 

XaRo 

This  is  a  general  relation  for  the  received  coherent  component  of 
pressure  for  any  surface  statistics,  degree  of  roughness,  and  obser¬ 
vation  point. 


Since  the  effective  extents  of  the  incident  illumination  on  the  surface 
are  much  larger  than  the  acoustic  wavelength,  P  takes  on  appreciable 
values  only  when  a  =  0,  b  =  0,  which  corresponds  to  the  specular 
direction.  (This  follows  upon  noting  that  P  in  (9)  is  a  double  Fourier 
transform  of  the  incident  pressure  'p.  .)  Therefore,  the  coherent 
component  is  appreciable  only  near  the  specular  direction. 


6 


For  a  perfectly  smooth  surface,  £  =  0,  and  the  first-order  CF 
equalsunity.  As  the  surface  roughness  increases,  the  magnitude  of  the 
CF  decreases,  thereby  causing  the  amplitude  of  the  coherent  component 
to  decrease.  The  coherent  component  is  not  identically  zero  for  a 
rough  surface;  however,  it  is  negligible  for  a  very  rough  surface. 

From  (10),  since  the  only  dependence  of  pc  (t)  on  time  is  via  the 
term  exp(i2.-fat),  the  spectrum  of  the  coherent  component  must  be 
an  impulse  at  frequency  fa  ,  just  like  the  transmitted  signal  spectrum. 
However,  the  amplitude  and  phase  shift  of  the  coherent  component  de¬ 
pend  on  the  physical  locations  of  the  transmission  and  observation 
points  and  the  degree  of  surface  roughness. 

The  scatter  component  of  the  received  pressure  is  defined  as  the 
remainder 

ps(t)  =p(t) -pr(t)  =  p(t) -p(t)  ,  (11) 

which  has  zero  mean.  In  order  to  evaluate  the  mean-square  value?  of 
the  scatter  component,  we  note  that 

|Ps(t)|2  =  ip(t)|2-|pc(t)|2=  |p(t)|2-  |Fw|2  .  (12) 


That  is,  the  mean  square  value  of  the  scatter  component  ps  is  equal 
to  the  mean  square  value  of  the  total  received  pressure  p  less  the 
squared  magnitude  of  the  coherent  component  p^  .  Using  (3).  we  obtain 


/Ojfdij dy j  dx2  dy2  p.(x,,y 


,)?*(x2,y2)  • 


(13) 


exp[ikaa(xj-x2)  +ik>b(y,-y2)lf^(kac,-kac;x,-x2,yJ-y2,0)  , 


where 


f^(kac, -kac;u,  v,  t  )  =  exp[ikac£(x,  y,  t)-ikacc(x-u,  y-v,  t-r)| 

=  d£2  exp(ikac£,  -  ikacc,)q2(  c,.  <C,;u,  v,  r) 


is  the  second-order  CF  of  the  surface  wave  heights.  Since  the  surface 
is  homogeneous  and  stationary,  q2  is  the  second-order  PDF  of 
surface  wave  heights  and  depends  only  on  differences  in  position  and 
time.  If  we  let  u  =  x  -  x2,  v  =  yt  -  y2  in  (13),  there  follows 


!p(t)i2=| 


/ - \  /fdudv  J(u, 


v)  exp[ika(a«4bv)]f^(kac,-kac;  u, v,0)  ,  (15) 
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Me»n  magnitude-squared  value,  more  precisely. 
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where  J  is  the  autocorrelation  of  the  illumination  p.  incident  on  the 
surface: 

J(u,v)  =  ffdxdy  p.U.ylp.*  (x-u,y-v)  ,  (16) 

Before  forming  |ps(t)j2  from  (12),  we  first  note  from  (9)  that  we  can 
express 

|P(k#a,kab)|  2  =  ffdudv  J(u,v)  expCik^au+bv)]  (17) 


Therefore,  using  (12),  (15),  (10),  and  (17),  we  obtain 

lPs(t)|2  =  |P(t)|2  -  |PcU)|  2 


ffdudv  J(u,  v)  exp[iki(au  +  bv)l  •  (18) 

|f £ (k a c,  — k a  c;  u,  v,  0  )  —  |f^(kac)  |  2J 

This  is  a  general  relation  for  the  mean- square  value  of  the  received 
scatter  component  in  terms  of  the  autocorrelation  function  J  of  the 
incident  pressure  p.  on  the  reflecting  surface  and  on  the  first-  and 
second-order  CF's  of  the  surface-height  variations  £•  No  assumptions 
about  the  degree  of  surface  roughness  have  been  made. 


We  will  now  make  the  reasonable  assumption  that  the  effective  ex¬ 
tents  on  the  surface  of  the  incident  illumination  p  are  much  larger 
than  the  distances  at  which  the  surface  heights  are  statistically  dependent 
on  each  other.  Mathematically,  this  is  equivalent  to  assuming  that  the 
bracketed  difference  of  CF's  in  (18)  decays  to  zero  in  u  and  v  much 
sooner  than  J(u,v)  does.  (The  difference  in  CF's  goes  to  zero  as 
ju|,  |  v  J  -*  «  because  then 

f?(kac,-kac;u,  v,  t)  — ♦  |fc(kac)|  2)  (19) 

The  double  integral  in  (18)  then  is  virtually  unchanged  if  J(u,v)  is  re¬ 
placed  by  its  origin  value,  yielding 

|p  (t)  |  2  =  I - ^  J  (0,0)  ffdudv  exp  [ik  (au+bv)] 

V»R°'  (20) 

[f.(kac,-kac;u,  v,0)  -  |  f,(kac)|  2I  . 

This  double  Fourier  transform  on  space  variables  u  and  v  is  as 
far  as  the  analysis  can  be  carried  without  further  assumptions  on 
wave-height  statistics.  The  surface  statistics  needed  in  these  quantitative 
measures  of  the  coherent  and  scattered  components  are  contained  entirely 
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in  the  first-  and  second-order  surface  height  CF‘s.  Arbitrary  surface 
roughness  is  allowed. 


As  the  surface  roughness  decreases,  the  difference  in  CF's  in  (18) 
and  (19)  goes  to  zero  (because  C  •*  0).  In  this  case,  the  scattered  com¬ 
ponent  disappears,  as  indeed  it  must  for  a  smooth  surface. 


In  Appendix  C,  the  scattering  strength6  S  of  the  surface,  defined 
as  the  ratio 


average  scatter  ;ntensity  at  receiver  due  to  unit 
scattering  art  %  referred  to  unit  distance 

astray  .cident  intensity  on  surface 


(21) 


is  shown  to  be  given  by 


IP,(t)l2  Kj 

cQJ(0,0) 


(22) 


If  we  substitute  (20)  into  (22),  the  above  equation  becomes 

o  2  | 

S= —  —  //dudvexp[ilca(au  +  bv)]  • 

CQ  X.  *  (23) 

[f?(k,c,-kac;u,  v,0)  -  |  f  ^(k^c)  |  2]  . 

This  is  a  dimensionless  quantity.  The  dimensionless  factor  B2/cq 
depends  solely  on  the  geometry  of  the  experiment. 


For  a  surface  of  very  slight  roughness,  i.  e. , 

k  # c  max  |  t(x,  y,  t) !  <<  1  , 

we  can  approximate  the  exponential  in  (3)  as 

expUktc  t(x,y,  t)]  =1  +  ikac  £(x,y,  t)  . 

Then,  (3)  becomes 

ft(t)  =-i  exp[i(2nfat  -  kaRo>]  [P(kaa,  k^b) 

•  O 

+  i  k a  c //  dx  dy  f>.  (x,  y)  exp li  k a  (ax  +by)l  £(x,  y ,  t)l  . 


(24) 

(25) 


(26) 


6The  term  "scattering  strength"  is  used  in  this  report  as  an  intensity  ratio  and  is  not  convened 
to  decibels. 
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The  first  term  of  (26)  is  the  coherent  component,  and  the  second  term 
of  (26)  is  the  scatter  component.  This  equation  for  the  coherent  com¬ 
ponent  is  a  special  case  of  (10),  where  CF  fj(k,c)  is  approximately 
unity.  Thus,  the  general  results  (10)  and  (18)  can  be  reduced  to  special 
cases,  including  very  slight  roughness  (or  very  rough  surfaces),  as 
desired,  by  appropriate  choice  of  CF's.  Also  (23)  is  a  general  relation 
for  the  scattering  strength,  which  is  applicable  to  any  degree  of  rough¬ 
ness  and  surface  statistics. 

3.  REFLECTION  FROM  A  NARROW-BAND,  TIME-VARYING  SURFACE 

In  this  section,  we  restrict  consideration  to  the  case  where  the 
surface-height  variation  at  each  point  of  space  is  a  narrow-band  function 
of  time  [10,  pp.  347-348  and  Section  8.5];  this  is  the  case,  for  example, 
when  the  sea  surface  is  characterized  as  swell  jlO,  Section  1.2]  . 

3.  1  NARROW-BAND  COMPONENTS  OF  RECEIVED  PRESSURE 
WAVEFORM 

The  surface  height  for  this  case  of  narrow-band  variation  can  be 
represented  as 

<(x,y,t)  =h(x,y,t)  cos[2fffst  +  0U,y,t)l  ,  (27) 

where,  for  fixed  position  x,  y,  the  amplitude  h  and  phase  0  vary 
slowly  with  time  t  in  comparison  with  cos(2irfst).  The  center  fre¬ 
quency  of  the  surface-height  variation  is  fs  .  If  we  substitute  (27)  into 
(3)  for  the  received  pressure  and  use  the  expansion  [ll,  8.5114] 

exp[ikach  cos(2rrf +  0)1 

OO 

=  2  i"  J  (k  ch)  explim(2?rf  t  +  0)1  , 

in  oo 

the  received  pressure  can  be  represented  as 

oo 

p(t)=  2  A  (t)  exp[i2>r(f  +mf  )  tl  , 

r  m=-®  “  as’ 

where 

Am(t)  =  -i  exp(— ik^) -JL  i“j;fdxdypi(x,y)  J  Jk.chU.y,  t)]  • 

A,Ko  W°) 

e*plika(ax+by)  +  im  0(*,y,  t)l  . 


(28) 

(29) 
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Since  h  and  9  vary  slowly  with  time,  so  also  does  Am(t).  However, 
if  k  ch  »  1  or  m  »  1,  the  rate  of  variation  of  A  (t)  is  much 
faster  than  that  of  h  or  9.  This  may  be  seen  by  noting  that  Bessel 
function  J  (x)  is  an  oscillatory  function  of  x;  thus,  Jm[kachj  can  go 
through  many  cycles  of  variation  while  h  goes  through  but  one,  if 
k  ch  »  1.  Also,  if  m  »  1,  exp(im9)  varies  much  faster  than  9. 

a  max. 

The  upshot  is  that  Am(t)  of  (30)  will  contain  a  considerably  higher 
frequency  content  than  either  h  or  9  if  either  of  the  above  conditions 
are  satisfied. 


Equation  (29)  expresses  the  received  pressure  waveform  as  a  sum 
of  narrow-band  components  (if  k  ch  and  m  are  not  much  larger 
than  unity),  with  the  mth  component  being  centered  at  frequency  f,  +  mfs 
That  is,  the  received  pressure  has  spectral  loLes  displaced  from  the 
transmitted  acoustic  frequency  fa  by  multiples  of  the  surface  center 
frequency  fs  ,  In  addition,  there  is  the  coherent  component  at  fa  .  This 
behavior  is  depicted  in  Fig.  2. 


-COHERENT  COMPONENT 


ZZ..JK 


fa+f> 


Fig.  2  -  Received  Acoustic  Specuum 


Asa  special  case  of  (30),  consider  a  traveling  sinusoidal  surface 
of  fixed  amplitude: 

h  (x,  y,  t)  =  h  , 

1  <31> 
(9(x,y,t)=c1x  +  c2y  +  c3  - 

Cj  and  c2  are  related  to  the  direction  of  travel  of  the  surface  wave 
and  its  velocity  of  propagation,  and  c?  is  the  phase  at  x  =  y  =  0.  Then, 
using  (9),  (30)  becomes 

g 

Ara(t)  =-i  exp(-ikaR0)^-~.  i“Jm(kach)  eq>(imc3)P(kaa  +  mcI,kab  +  mc2)  ,  (32) 

•  O 

which  is  independent  of  time.  This  case  has  been  investigated  previously 
[?,  Appendix]. 
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Now  if  the  surface  is  actually  random,  the  quantities  h  and  c3, 
above,  could  be  considered  as  random  variables,  with  the  PDF  for  c? 
uniform  over  a  2it  interval.  In  this  case,  the  mean  ( ensemble  average) 
of  the  mth  complex  amplitude  Am(t)  in  (32)  is  zero  for  m  /  0;  thus, 
there  are  no  coherent  components  in  any  of  the  sidebands.  For  m  =  0, 
an  assumption  on  the  PDF  of  h  is  necessary  to  evaluate  the  ensemble 
average  Jm(kach).  Forexample,  for  a  Rayleigh  PDF  of  surface  heights , 

p-cxp(“p)‘h>0  •  (331 

where  a  is  the  rms  wave  height,  the  mean  value  of  A0  (t)  [ll,  6.631  4] 
is  _ 

A  (t)  - -i  exp  (-1  k  R  )  - P  (k  a,  k  b)  exp  (-'i  k  2c2o2)  .. 

°  K\  '  3  (34) 

However,  (34)  coupled  with  (2^)  is  seen  to  be  but  a  special  case  of  (10) 
when  £  is  first-order  Gaussian: 

f  ^(k a  c)  -  exp(-‘4k2c2a2)  .  (35) 

Thus,  the  approach  given  in  Subsection  2.2  is  a  very  powerful  one  for 
evaluating  the  coherent  component,  and  includes  numerous  special 
cases.  Nevertheless,  (29)  and  (30)  are  useful  for  lending  insight  into 
the  spectral  behavior  cf  the  received  pressure  waveform  for  a  narrow- 
band,  surface -height  variation.  They  could  be  used  as  che  starting  point 
for  the  theory  to  be  developed  in  the  next  subsection,  but  a  more  com¬ 
pact  approach  has  been  utilized  there. 

3.2  SPECTRA  AND  POWER  OF  INDIVIDUAL  NARROW-BAND 
COMPONENTS 

The  received  pressure  waveform  p(t)  was  given  in  (3)  and  the 
coherent  component  (the  mean  of  the  received  pressure)  in  (10).  We 
now  wish  to  evaluate  the  correlation  function  of  p,  and  then  the 
spectrum,  for  a  narrow-band,  surface-height  variation.  The  correlation 
function  Rp  of  p  is,  by  a  simple  generalization  of  (13)  through  (15), 


Rp  (A  =  p(t!  p*  (t-  r) 

exp(i2zrfaf)  ff  dudvj(u,v)  exp[ik#(au +bv)lf  ^(kac,-kac;u,  v,  A 


7 .Since  the  surface  process  and,  therefore,  the  received  acoustic  process  .ire  noi  ergodic  in  this 
case,  this  average  is  not  equal  to  the  time  average  over  individual  member  functions. 
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(36) 


p(t)  p' 


erp(i2n'fa  r)  ff  dudv  J  (u,  v)  exp[ika(au  +bv)l 


[f?(kac,  -kac;U, v,  r)  -  )  f  ^(k a c)  |  2]  , 

where  we  have  used  (10)  and  (17).  We  define  the  covariance  function  of 
p  as  the  correlation  function  of  the  ac  component  (i.  e.  ,  scatter  com¬ 
ponent)  of  p: 

R,M  =  PSM  ps*  U-rt 

=  (p(t)  —  p  (t)]  [p*  (c  —  r)  -  p*  (t  -  r)] 

=  e*P  (*2»f#»')  f/du  dv  J  (u,  v)  exp[ika(au  +bv)]  •  ($7) 

(f5(kae, -kae;u,v,  r)  -  |f?(kac)|  2]  . 

Equation  (37),  a  general  relation  for  the  correlation  function  of  the 
scatter  component,  reduces  to  the  mean-square  pressure  of  (18)  for 
r  =  0. 

Again  we  make  the  assumption  that  the  effective  extents  on  the 
surface  of  the  incident  illumination  p.  are  much  larger  than  the  dis¬ 
tances  at  which  the  surface  heights  are  statistically  dependent  on  each 
other,  andget  the  approximation  (See  (18)  through  (20))  for  the  correlation 
of  the  scatter  component: 

Rs(r)  =  {■  J  exp(i2wfar)  J  (0,0)  ff  oudveiplik  (au+bv))  w 

\A»Ro/  *  (38) 

[fj(kac, -kae;u.  v,  r)  -  |f^(kac)|2l  .. 

At  this  point,  we  make  an  assumption  about  the  statistics  of  the  surface 
heights,  namely,  that  the  second-order  PDF  of  surface  heights  is  joint 
Gaussian  [lO,  pp.  343-345].  Then  [l2,  Eq.  (8-23)], 

f  £  (^a c>  ~k,c:u» v*  *)  =  expl— k2 c2<*2(l  —  p(u,v,r)}l  ,  (39) 

where  o  is  the  rms  wave  height  and  p  is  the  normalized  spatial- 
temporal  correlation  function  of  surface  heights ,  assumed  homogeneous 
and  stationary: 

v2p(u,  v,  r)  =  <f(x,y,  t;  £(x-u.y-v,  t-r)  ..  (  ">) 


If  we  substitute  (39)  into  (38)  and  use  (19),  there  follows  for  the 
correlation  function  of  the  scatter  component 

R  (r)  -  (— exp(i2>rf  r)  J  (0,0)  exp(-/32)  |fdudvexp[ik  (au+bv)]  • 

\AaRc/  (41) 

[expl$2p(u,  v,  r)| -l]  , 

where  we  have  defined  the  surface  roughness  parameter  (2,  p.  82,  Eq.  (10)] 

0=kac  a  ,  (42) 

The  mean-square  value  of  the  scatter  component  for  the  Gaussian  surface 
is  given  by  substituting  r  =  0  into  (41);  then  only  knowledge  about 
p(u,v,  or.  the  normalized  correlation  between  two  separated  surface 
heights  at  the  same  instant  of  time,  is  required. 

Thus  far,  the  surface  correlation  function  p(u,v,  r)  has  been 
general.  The  spectrum  Gs  of  the  scatter  component  is  given  by  the 
Fourier  transform  of  (41), 

Gg(f)  = /drexp(-i2fffr)  R s ( r)  ,  (43) 

and  generally  can  be  numerically  evaluated  by  a  Fast  Fourier  Transform, 
(FFT).  Thus,  the  received  scatter-pressure  spectrum  could  be 
evaluated  from  (41)  and  (43)  for  any  degree  of  surface  roughness,  with 
no  assumptions  about  narrow -band  surface  variation  and  narrow-band 
components  of  the  received  pressure.  However,  a  double  integral  and 
an  FFT  is  involved. 

The  approach  taken  in  this  report  is  to  specialize  to  the  case  of  a 
narrow-band,  surface-height  spectrum.  Then,  p  takes  the  form  (See 
Appendix  D  ) 

p(u,  v,  t)  *  Re  i  g(u,  v,  r)  exp(i27rfsr)j 

=  i  g(u,  v,  t)\  cos  [2irfsr  +  arg  ig(u,  v,  r)|]  , 

where  jj(u,v,  r)  is  the  complex  envelope  of  p(u,v,  r),  and  varies 
slowly  with  r  as  compared  to  exp(i2tr  fsr). 

Upon  substitution  of  (44)  into  (41),  and  using  the  expansion  [ll,  8.511  4 
and  8.  406  3] 

OO 

exp(xcosci)-  S  Im(x)  exp(imd>)  ,  (45) 
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we  get 


R  At)  =  /  ..R  A  e*p(i2fff  r)  J(0,0)  e*p(-j82)  £  //dudv  (46) 

V‘.Ro/  m=““ 

exp  [i  k  a  (au  +  bv)]  |jm  (/32  |  £(u,  v,  r)  j )  -  SoJ  exp(i27rmfsr  +  imarg|£(u,  v,  r)l)  , 
where  S  is  the  Kronecker  delta. 

OID 

The  spectrum  of  the  scatter  component  is  obtained  by  Fourier 
transforming  (46): 

G*(f)=(ir)2  ,(0,0)«io.T»(f"f«_mf*)  *  (47) 

where 

Tn  (f)  = /dr  exp(— i2»rfr)  exp  (-/S2)  /f  dudv  exp  (ika(au+bv)J  •: 

(48) 

Thus,  the  scatter  spectrum  is  composed  of  spectral  lobes,  or  sidebands, 
centered  at  frequencies  equal  to  the  acoustic  frequency  f,  plus 
multiples  of  the  surface  center  frequency  f$  ,  as  anticipated  by  (29) 
and  (30)  and  shown  in  Fig.  2.  The  mth  order  sideband  is  defined  as 
that  spectral  lobe  centered  at  frequency  fa  +  mf  s .  The  zeroth  order 
sideband  is  centered  at  f,  ,  but  is  spread  in  frequency;  it  is  distinct 
from  the  coherent  component  of  the  received  pressure,  which  has  a 
delta  function  at  frequency  f, . 

For  a  given  surface  spatial-temporal  correlation  function  p, 
observation  point  (a,  b,  c),  and  surface  roughness  p,  (48)  can  be 
numerically  evaluated  (by  a  double  integral  and  FFT)  for  the  spectrum 
of  the  mth  sideband  component  of  the  received  scatter  pressure.  The 
total  power  in  the  mth  sideband  is  obtained  by  integrating  the  mth  term 
in  (47)  over  all  frequencies,  and  is  denoted  by 


jlm(/32|g(u,v,r)|)  -6omj  exp(imargle(u,  v,  r)J) 


i  i  /  n  \  2 

|p1B(t)|2  -  (-■■■■■  )  J  (0,0)  exp (— /S2)  //dudv exp [ika(aa+bv)l 

V  *^°/ 

e(n,  ▼,  0)  | )  -50J  exp(imarg{g(u,  v,0)l)  . 


(49) 
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The  scattering  strength  of  the  mth  sideband  is  defined  in  a  manner 
similar  to  (C - 12)  in  Appendix  C  as 

I  p  }  2  r  2  g2  j 

S  °  .  0  =  —  — «•  e*p(-/32)//dudve«p[ilt.(*u+bv)]  • 

CqJW.O)  cq  \  (50) 

[lB  (fi2  I  V,  0)  I )  -  *.J  e*p  (i  m  arg  Igju,  v,  0)0  , 

The  scattering  strength  Sa  depends  on  the  spatial -temporal  correlation 
function  p.(u,v,  r)  at  zero  time  delay  (  r  =  0). 


3.  3  SPECIAL  FORM  OF  SURFACE  CORRELATION 

Thus  far,  the  spatial-temporal  surface  correlation  function  p(v ,  v,  r) 
has  been  general,  except  for  the  narrow-band  assumption.  We  now  wish 
to  specialize  to  a  particular  form.  Note  first  that  if  z  is  purely  real, 
but  perhaps  negative, 

Iffl(  |z|)  exp(imarglzl)  =IB(z)  .  (51) 

Now  consider  that  correlation  function  p  has  the  form 

-  <52> 

That  is,  for  a  fixed  delay  r  ,  contours  of  iso-correlation  values  are 
elliptical.8  (There  is  no  need  to  consider  a  rotated  ellipse  if  the  x,  y 
axes  are  aligned  with  local  surface  directional  properties. )  Distances 
Lx  and  Ly  are  the  (correlation)  distances  in  the  x  and  y  directions, 
respectively,  at  which  the  correlation  is  down  to  a  specified  fraction 
(e.  g. ,  l/e)  of  its  peak  value. 

Several  important  special  cases  of  surface  correlation  can  be 
subsumed  by  (52).  For  example,  if  the  correlation  distances  Lx  and 
Ly  are  equal,  the  surface  is  isotropic.  However,  if  one  of  the  corre¬ 
lation  distance  is  infinite,  the  surface  correlation  function  is  one¬ 
dimensional;  i.  e. ,  it  depends  on  only  one  of  the  variables  u,v.  Further 
specialization  of  the  "one-dimensional"  surface-correlation  function 
would  be  a  periodic  "one-dimensional"  surface -correlation  function, 
and  as  a  particular  case  of  the  latter,  a  sinusoidal  surface-correlation 


8  See  for  example,  Ref.  13,  p.  81,  where  experimental  results  of  this  form  have  been  obtained. 
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function.  Thus  (52)  is  particularly  useful  for  detailed  investigation 
because  it  retains  a  great  deal  of  generality  and  is  a  reasonable  approx¬ 
imation  to  many  practical  situations. 


Using  (D-2)  and  (52),  we  note  that 

p(u,  v,  -t)  =  p(u,  v,  t)  . 


(53) 


That  is,  surface  correlation  p  is  even  in  r  for  any  u,v,  under 
assumption  (52).  Employing  (D-15)  through  (D-17),  wehavefor  this  case 
of  (52), 


(54) 


If  we  combine  (50),  (51),  and{54),  the  mth  order  scattering  strength  is 


n2  t 

- - y  e*p(-/S2)  Jjfdu  dv  exp[ika(au  +bv)] 

CQ  A. 


(55) 


Since  I_Jx)  =  Im(x),  the  scattering  strengths  of  the  pair  of  symmetrically 
located  sidebands  at  fa  -  mf ,  and  fa  +  mf,  are  equal,  under  the 
assumption  (52);  this  is  not  true  for  general  p.  When  we  let 

u  =Li  t  cos  0  , 


v  =  L  t  sin  0  , 
y 


(56) 


in  (55)  and  perform  the  integration  on  6,  the  scattering  strength  becomes 


sn=~^2-«p(-/32)7  *‘U«>[lB(|8V'>>-5oJ  • 


Q  Xl 


where 


[HkiV«2L2+b2L 


Recollect  from  (42)  that 


^=k  Cff 


(57) 

(58) 

(59a) 


The  fundamental  parameters  a  and  p  are  dimensionless  and  are  basic 
measures  of  surface  behavior  in  the  horizontal  and  vertical  directions, 
respectively.  Tne  former  is  related  to  the  correlation  distance  of  the 
surface,  whereas  the  latter  depends  on  the  surface  roughness. 


17 


The  quantity  a  in  (58)  can  take  values  in  the  range  from  zero  to 
very  large  values  compared  with  unity.  For  example,  in  the  specular 
direction,  a  =  b  =  0  and  a  =  0.  However,  for  other  directions,  if  the 
ratio  of  correlation  distance  to  acoustic  wavelength  is  much  larger 
than  unity,  a  takes  on  large  values.  Of  course  a  can  not  become 
arbitrarily  large  by  letting  the  correlation  distances  increase,  because 
they  have  been  assumed  less  than  the  illuminated  extents.  (See  text 
following  (18).)  The  quantity  (3  in  (59)  can  take  values  in  the  range 
from  zero  to  values  of  the  order  of  10  without  violating  the  conditions 
stated  in  the  paragraph  following  (30).  Thus,  for  a  smooth  surface, 
P  =  0,  whereas,  for  a  rough  surface,  p  can  take  on  very  large  values. 
But,  in  this  latter  case,  the  sideband  components  in  Fig.  2  of  the  received 
scatter  pressure  would  spread  out  significantly  infrequency  and  over¬ 
lap  each  other.  These  sidebands  can  not  be  separated  at  the  receiver 
by  a  filter,  and  the  scattering  strength  Sa  would  lose  its  meaning. 

It  would  be  necessary  to  resort  to  the  general  case  given  in  (41)  and 
(43)  for  the  received  scatter  spectrum  to  find  how  much  power  lies  in 
a  particular  spectral  band. 

The  scatter  correlation  function  for  the  elliptical  spatial  correlation 
case  is  obtained  by  substituting  (52)  into  (41): 

R  (r)  =  2»L  L  (-1^  e*p(i2»rf  r)  J(0,0)  exp(-/32)  • 

\W  (59b; 

/  dr  r  Jo  (at)  [eip  {/32  pl  (r,  r)l  -  l] 

O 

However,  for  our  purposes  in  this  report,  in  order  to  retain  physical 
significance  and  interpretation  for  scattering  strength  Sffl  ,  we  consider 
a  and  p  upper-limited  to  values  of  the  order  of  10. 

For  the  specular  direction,  (57)  becomes  (See  Appendix  A) 

S»*=CQ(k.L«)(k.L,)U-^)  ’  <60> 

where 

ua  (@)  =  (^r1  txp(-^J)  /  dt  r  [lB(^P2W)  -  SOB]  .  (61) 

The  function  Ua  (P)  in  (61)  depends  on  surface  roughness  p  and 
the  form  of  the  spatial  correlation  function  p2  .  In  Section  5,  plots  of 
U  (p)  versus  P  for  several  forms  of  spatial  correlation  p2  are 
given,  including  exponential  and  Gaussian  spatial  correlation  and 
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exponential  and  Gaussian  modulation  of  a  sinusoidor  Bessel  function  for 
the  spatial  correlation  The  factor  preceding  Um(P)  in  (60)  is  held  con¬ 
stant  for  these  plots;  thus,  the  geometry,  acoustic  frequency,  and  corre¬ 
lation  distances  are  considered  fixed.  Each  plot,  therefore,  measures 
the  dependence  of  the  mth  scattering  strength  on  surface  roughness 

P  as  the  rms  wave  height  a  is  varied. 

For  directions  other  than  specular,  a  /  0,  and  a  change  of  variable 
in  (57)  yields 


S 


m 


V„(a,/3) 


(62) 


where 

v.(a,P)  exp(-£2)  /  ds  s  J0 (s)  [lm  (02  p2 (s/a))  -  Som]  . 

The  function  Vm(a,P)  in  (63)  depends  on  a,  p,  and  the  form  of  the 
spatial  correlation  function  p2  .  If  the  factor  preceding  V.,  in  (62)  is 
held  constant,  the  geometry  and  the  ratio  of  correlation  distances  must 
be  considered  fixed.  Then,  the  function  VB  measures  the  dependence  of 
scattering  strength  on  a  as  the  correlation  distances  La  and  L.y  are 
varied  (although  their  ratio  is  fixed)  and  on  P  as  the  rms  wave  height  a 
is  varied.  Plots  of  Vm  (a,  P)  versus  a  and  p  for  several  forms  of 
p2  are  also  presented  in  Section  5. 


For  the  particular  surface  correlation  function  form  assumed  in 
(52),  the  general  expression  for  the  correlation  function  of  the  received 
scatter  pressure  waveform  given  in  (41)  takes  a  special  form.  It  is, 
using  (56), 


R,(r)  =  B2(k,Lj[)  (k. 


Ly> 


J  (0,0) 


exp(i2n-far)  exp(-/32)  /  dr  r 

O 


Jo  (ar)[esp  I/32  Pj  (r,  r)l  -  l]  . 


(64) 


Thus,  only  a  single  integral  need  be  evaluated  in  order  to  obtain  the 
scatter  correlation  function.  The  scatter  spectrum  follows  from  (43). 
Arbitrary  surface  roughness  is  allowed.  The  scattering  strength  for  the 
total  scatter  power  is,  using  (C-12)  and  (54), 
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(65) 


*,«»  Rp 
cQJ  (0,0) 


B2  no 

=  —  (k#Li)  (k^L  )(2ir)~l  e*p  (-/32)  /  dr  r  • 

C  Q  ”  O 


J  o  (a  r)  [  e«p  1/3  2  p  2  (r)  I  -  l]  . 

Again,  a  single  integral  must  be  evaluated  for  the  scattering  strength, 
and  arbitrary  surface  roughness  is  allowed. 


.1.4  NARROW-BAND  SIDEBAND  SPECTRA  FOR  SMALL  ROUGHNESS 


The  spectrum  of  the  received  scatter  pressure  was  represented 
as  a  sum  of  narrow-band  lobes  in  (47)  and  (48)  for  general  surface 
roughness  p  and  narrow-band  correlation  g.  (See  Fig.  2. )  For  small 
P,  a  useful  approximation  to  (48)  can  be  made:  first  note  [l  1 ,  8.445] 


that 


(66) 


When  we  employ  this  approximation  in  (48)  and  use  (D-7),  there  follows 
To  (f)  /94  exp(-/32)  //dudvexp[ik((au  +  bv)]  /dw  ^(u,  v,  w)^*  (u,  v,  w-f)  , 

T  t  (£)  =—  /32  e*p  (~/32)  /fdu  dv  e*p[i  kB(au +bv)]  ^(u,  v,  f)  , 

T2(f)  /94  e*p(-/32)  //du  dv  exp[i  k  a  (au  +  bv)]  /  dw  ^{u,  v,  w)  ^(u,  t,  f-w)  . 
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The  largest  term  is  Tt  (f),  which  is  proportional  to  p2  and  to  the 
double  Fourier  transform  of  the  low -pass  spectrum  £  of  the  surface 
correlation.  Tj  (f)  is  not  simply  proportional  to  g(0,0,f),  the  low  - 
pass  spectrum  at  a  point  of  the  surface -height  variation.  Tq  (f)  is  pro¬ 
portional  to  P4  and  is  related  to  the  autocorrelation  of  the  low -pass 
spectrum  g.  T^  (f)  is  proportional  to  p4  and  related  to  the  convolution 
of  the  low-pass  spectrum  g. 


For  the  special  case  when  the  spatial-temporal  surface  correlation 
function  is  separable  in  space  and  time  variables. 
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(68) 


P  (a*  V,  r)  =  p}  (u,  v)  p4  (r) 

g(u,v.f)  is  also  separable, 

g(u,  v,  f)  =  p3  (u,  v)  g4  (£)  ,  (69) 

and  the  frequency-dependent  terms  of  (67)  are  given  by 

/dwj54(w)  g4  (w-f)  ,  (70a) 

g4(f)  ,  (70b) 

/ dwg4(w)  g4(f-w)  ,  (70c) 

respectively.  Since,  from  (69), 

g(0,0,f)  =  g4(f)  ,  (71) 

g4(f)  is  proportional  to  the  surface-height  spectrum  at  a  point  and, 
therefore,  is  real.  Equation  (70)  states  that  the  zeroth-order  scatter 
sideband  spec trum  Tq  (f)  is  directly  proportional  to  the  autocorrelation 
of  the  wave-height  spectrum  and  is  even  about  f  =  0  (which  corresponds 
to  f  =  f,  in  the  received  acoustic  spectrum).  Also,  the  first-order 
sideband  spectrum  is  directly  proportional  to  the  wave-height  spectrum, 
and  the  second-order  sideband  spectrum  is  proportional  to  the  convo¬ 
lution  of  the  wave-height  spectrum  and  need  not  be  even  about  f  -  0 
(which  corresponds  to  f  =  fa  +  Zig  in  the  received  acoustic  spectrum). 

The  mth  order  sideband  Tn  (f)  is  proportional  to  p2”  for  small 
p  and  |m|  >  1.  This  sideband  involves  higher  order  convolutions  of 
the  low-pass  spectrum,  thereby  causing  significant  spreading  in  fre¬ 
quency.  This  is  consistent  with  the  observations  made  in  the  text 
following  (30). 


3.  5  RMS  BANDWIDTH  OF  SIDEBANDS  FOR  SMALL  ROUGHNESS 

It  is  of  interest  to  know  quantitatively  the  amount  of  frequency 
spreading  that  each  sideband  undergoes.  One  simple  measure  of  this 
spreading,  short  of  evaluating  the  actual  sideband  spectrum,  is  the 
rms  bandwidth  of  each  sideband.  We  start  by  defining  the  rms  band¬ 
width  Bs  of  the  surface-height  spectrum  g  for  the  special  case  of 
(68)  according  to 
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(72) 


b^-L^P.’0,0  -HULbM.fdtt*  -(I)  , 

'  /dfglO.O.f)  / df  g4  (f) 

where  we  have  used  Appendix  D  and  the  fact  that 

P,( 0,0)  =p4( O'  -  1, 

since  p(0,  0,  0)  =  1. 


(73) 


For  small  roughness,  the  spectra  of  the  zeroth  and  first- and  second- 
order  sidebands  are  given  in  (70).  From  (70b),  it  follows  immediately, 
using  (72),  that  the  rms  bandwidth  Bj  of  the  first-order  sideband  is 
equal  to  Bs ,  the  rms  surface-spectrum  bandwidth.  The  zeroth  order 
rms  bandwidth  Bo  is  available  from  (70a)  as 


/ dff2  fd w  g  (w)  g  (w-f) 

B2  - - — - — - =  2B2 

f  df  /  dw  g4(w)  g4  (w-f) 

utilizing  the  fact  that 

/dff  g4 (f)  =  0  , 


(74) 


(75) 


from  (D-5),  (D-6),  and  (69).  Similarly,  it  may  be  shown,  in  general, 
that  the  rms  bandwidth  of  the  mth  order  sideband  is  given  by 


B 


ffi 


m 


ID  ^  0 


(76) 


if  0  «  1.  Thus,  the  zeroth  and  second-order  sidebands  are  V?  wider 
than  the  surface  spectrum,  and  the  higher  order  bandwidths  increase 
as  the  square  root  of  the  order  number.  For  the  small  roughness  case, 
the  correlation  distances  or  direction  cosines  do  not  enter  into  this 
relation. 


4.  RELATION  OF  SURFACE  CORRELATION  FUNCTION  TO 
DIRECTIONAL  WAVE  SPECTRUM 

The  results  in  the  previous  section  require  specific  forms  for  the 
spatial-temporal  surface  cor  relation  p  for  their  numerical  evaluation. 
In  this  section,  we  shall  use  the  relation  between  the  surface  correlation 
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function  p  and  the  directional  wave  spectrum  for  a  wind-generated 
sea  composed  entirely  of  gravity  waves,  and  specialize  to  several 
cases  of  interest. 


4.  1  GENERAL  RELATIONS 


For  a  homogeneous  surface,  the  spatial-temporal  surface  corre¬ 
lation  function  p(u,v,r)  can  be  expressed  in  terms  of  the  directional 
wave  spectrum  A2((jl,^)  [l3,  Part  8;  10,  Chapter  8j- 

a 2  p(u,  v,  r)  =  ff  dp  dv  %2  (p,  v)  cos [up  +  vu  —  g^  (p2  + 1/2)^  r]  ,  (77) 

where  g  is  the  acceleration  of  gravity  and  p  and  u  are  the  wave 
numbers  in  rectangular  coordinates.  In  particular,  the  temporal  corre¬ 
lation  of  the  surface-height  variation  at  a  point  is  given  by 

<?2  p(0,  0,  r)  =  ff  dp  dv  ft2  (p,  v)  cos\gVl  (p2  +  v2)'4  t]  ..  (78) 


An  alternate  representation  of  the  directional  wave  spectrum  in  terms 
of  the  polar  form  A2  .p,v)  is  often  used: 

g'^  A2  (g,/*  (p2  +  i'2)'4,  tan-1  v/p) 


A2  (p,  i')  =• 


(79) 


2  (p2  +  v2)y> 

If  we  substitute  (79)  into  (78)  and  make  the  substitutions  p  =  (2irf)2  cos  9/g, 
v  -  (2trf)  sin9/g,  the  surface-height  temporal  correlation  becomes 


a2  p(0,0,  t)  =  2rr  f  dfcos(27rfr)  /  dd  A2  (2rrf,  6) 


(80) 


The  (double-sided)  surface-height  spectrum  (at  a  point)  is  denoted  by 
4>(f)  and  is  the  Fourier  transform  of  (80): 

n 

<t>(£)  =  / dr  exp  (-i2?rfr)  o2  p(0,  0,  r)  =  n  f  dd  A2 (277-  [  f  |,0) 


=  (2rr>4  g’2  |f \>fAaV(**ltleo*e  ,  4*2f2sinfl\ 

V  8  g  / 


(81) 


The  last  form  follows  from  (79). 


4.2  ELLIPTICAL  SURFACE  CORRELATION 

The  results  above  are  for  general  directional  wave  spectra.  Here 
we  will  consider  a  special  case  of  particular  interest.  We  assume  first 
that  the  directional  wave  spectrum  ^2  possesses  180°  symmetry: 
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A  2  (— fi,  -u)  =  X2((i,v) 


(82) 


This  is  done  for  mathematical  tractability,  rather  than  for  any  under¬ 
lying  physical  reason.  Then,  from  (?7),  p(u,v,  t)  is  even  in  r  for 
all  u,v: 

p(u,  v,  -r)  =  p(u,  v,  t)  . 

When  we  use  (D-15)  through  (D-17)  and  (77),  there  follows 

d (u,  v,  0)  =  p (u,  v,  0)  =  dp dv  o~ 2  A2  (p,  v)  cos  (ep  +  vu)  (84) 

But  this  ‘relation  can  be  inverted  by  a  double  Fourier  transform  to 
obtain  the  directional  wave  spectrum  £2(p,t/)  in  terms  of  p.(u,v,  0): 

a~2  A2(p,  t')  =  (2tt)'  2//dudvexp[-i(pu  +  i/v)]  p(u,v,0)  .  (85) 

We  have  used  the  symmetry  of  ^2  in  obtaining  (85). 

For  the  elliptical  correlation  form  assumed  in  (52),  we  had  in  (54) 


If  we  substitute  (86)  into  (85),  the  directional  wave  spectrum  takes  the 
form 


Thus,  the  directional  wave  spectrum  is  also  elliptical  if  the  surface 
correlation  is  elliptical  at  zero  delay.  (An  isotropic  surface  is  a  special 
case  of  (87).)  Equatic  i  (87)  may  be  expressed  compactly  as 

A 2  (17''  =  a2  L  L  (2ff)-1  f  drrp  (r)j  for)  ..  (88) 

z  *  y  0  1  ° 

For  any  particular  form  of  spatial  correlation  p2  ,  (88)  can  be  evalu  ¬ 
ated  numerically;  then,  (87)  yields  the  directional  wave  spectrum. 

Conversely,  if  the  directional  wave  spectrum  %2  is  elliptical,  so 
also  is  ft(u,v,0).  The  exact  relation  is  obtained  by  setting 
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A2(p,v)  =  A2  (y}(Lttl)U(Lyuy 


in  (84).  There  follows 


w 

l2.(u,  v,  0)  =  (a2  Lx  Ly)~ 1  2n  /  drj  rj  A 2(rj)  J£ 


or  more  compactly, 

p2(r)  =(ff2LsL  V*1  2»r/ djj  7  AjG;)  J0(n7)  .  (91) 

O 

Equations  (88)  and  (91)  are  a  Hankel-transform  pair  of  order  zero 
(j.4,  p.  136j .  Specification  of  either  p2  or  A2  determines  the  other. 

The  surface-height  spectrum  $(f)  for  the  elliptical  surface  corre¬ 
lation  is  obtained  by  substituting  (89)  into  (81): 


<Kf)  =  (2*r)4  ff 2  |f|3  /  d(9Aj  \  J(Lxcos0)2  +(Lrsiad)2 


S' 


For  the  special  case  of  an  isotropic  surface,  L,  =  Ly  =  L,  (92)  becomes 


<Mf)  =  (2*)5g'2  If |3  A2 

and  (91)  takes  the  more  familiar  form 


o2p2(t)  =  /  df[24>(f)l  J( 


¥r) 

■<r& 


(93a) 


(93b) 


Here,  2<t(f)  is  the  equivalent  single-sided  surface -height  spectrum. 


4.  3  EXAMPLES 


Four  examples  of  the  elliptical  surface  correlation  function  given  in 
(86)  will  be  considered.  The  first  is  exponential  decay  of  a  cosinusoidal 
spatial  variation: 
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p2( r)  =  exp(-Rr)  cos  (Qr),  r  >  0 


(94) 


For  Q  equal  to  zero,  simple  exponential  decay  of  the  spatial  corre¬ 
lation  is  realized.  R  is  of  the  order  of  unity.  When  vve  substitute  (94) 
into  (88),  it  follows  that  [ll,  6.623  2] 

{1+iQ/R  ) 

- (  (95) 

[(1+iQ/R)2  +  (17/R)2] 3/2) 

The  surface-height  spectrum  for  the  elliptical  surface  correlation  then 
requires  the  numerical  evaluation  of  (92).  For  an  isotropic  surface, 
I,x  =  Ly  =  L,  substitution  of  (95)  into  (93a)  yields  for  the  surface -height 
spectrum 


4>(f)  =  2no2 


(96) 


where  the  dimensionless  parameter  x  is  defined  as 

fgR  . 

and  the  dimensionless  function  $  is  defined  as 

1+iQ/R 


<DU)  =  |*!5Re< 


i_4 


(97) 


(98) 


(((1+iQ/R)2 

A  plot  of  the  frequency-dependent  term  $  is  given  in  Fig.  3  for  several 
values  of  Q/Rj<  1.  As  x  -*  0+,  the  Re  component  of  (98)  approaches 

1— (Q/R)2 

(99) 


[1  +(Q/R)2)2  ' 

Since  the  surface -height  spectrum  can  not  be  negative,  it  is  necessary 
(but  not  sufficient)  that  Q/R  <.1  in  order  for  (94)  to  be  a  validfo-m  of 
spatial  correlation.  The  spectrum  of  (96)  decays  as  | f l~ 3  for  large 
frequencies. 


Since  Q  was  upper-limited  to  a  value  of  unity  in  the  above  example, 
the  surface -height  spectrum  can  not  be  made  arbitrarily  narrow.  It  is, 
therefore,  of  interest  to  demonstrate  that  arbitrarily  narrow  spectra  are 
possible  through  a  slightly  modified  spatial  correlation,  namely,  ex¬ 
ponential  decay  of  a  Bessel  function: 
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Fig.  3  -  Surface-Height  Spectrum  for  Exponentially  Modulated  Cosine  Spatial  Correlation 
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p2  (r)  =  exp  (-Rr)  Jg  (Qt),  r  J>  0 


(100) 


This  function  also  oscillates  with  r,  but  decays  slightly  faster  than 
(94)  because  the  envelope  of  J0  (t)  decays  as  t“H  for  large  t.  If  we 
substitute  ( 100)  into  (88),  there  follows  [l 5,  pp.  314-316,  Eqs.  (11),  (14), 
and  (19)] 


A2(//)  =  a2LxLy  (2tt)~  1  [  drrexp(-Rx)  Jo (Qr)  J0  G;r) 


j  r  q-l  r  ,-i-H 

l+(±-l)2  e (•)  , 

n1  R2  |_  \R  Rj  J  |_  VR  R/  J 


(101) 


where  E(*)  is  a  complete  elliptic  integral  of  the  second  kind  with 
modulus  [l6,  Chapter  17;  see  especially  p.  590] 


or  parameter 


(102) 

(103) 


The  surface -height  spectrum  follows  upon  substitution  of  (101)  into 
(92)  for  the  elliptical  surface  correlation  case.  Since  Aj  (q  )  of  (101) 
is  positive  for  all  choices  of  Q  and  R  [16,  p.  609;  m  of  (103)  is 
always  less  than  unity],  the  surface-height  spectrum  4»(f)  of  (92)  is 
nonnegative;  therefore,  we  consider  the  spatial  correlation  form  in 
(100)  for  arbitrary  values  of  Q  and  R. 


Rather  than  evaluate  (92)  numerically,  we  restrict  attention  here 
to  the  isotropic  surface  and  use  (93a)  to  obtain  the  surface -height 
spectrum: 


where 


and 


4>(f)  =  2no2  J—  $(x)  , 
IgR 


x  =  2nf 


(104) 


(105) 


sw.i.ixpfi^x2-^)2]'1^^  E(.) .  <io6> 
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The  parameter  m  of  the  complete  elliptic  integral  is  given  by 

i-i 


ent  function  $(x)  is  p! 


(107) 


The  frequency-dependent  function  4>(x)  is  plotted  versus  x  in  Fig.  4 
for  several  values  of  Q/R.  The  spectrum  decays  as  |  f  |~ 3  for  large 
frequencies.  The  narrow-band  character  of  the  surface-height  spectrum 
for  large  Q/R  is  evident. 


The  third  example  of  spatial  surface  correlation  is  Gaussian  decay 
of  a  cosine9  : 


p2(f)  =  exp(-R2t2)  cos(Qr),  r_>0 


(108) 


Substituting  (108)  into  (88),  we  obtain 


.•(£•)  ■  U09) 

The  surface-height  spectrum  for  an  isotropic  surface  follows  upon 
substituting  (109)  into  (93a): 


where 


4>(f)  =  2na  2 


(110) 


dll) 


and 


<&(x)  =  jx|3/  dt t exp (-t2)  cos J  (x2t)  .  (112) 

This  function  is  plotted  in  Fig.  5  for  several  values  of  Q/R.  (The  method 
for  evaluating  (112)  is  given  in  Appendix  E.  For  Q/R>  1.  848,  $(x) 
goes  negative  at  the  origin,  thereby  invalidating  (108)  in  that  range.) 


The  fourth  example  of  surface  correlation  to  be  considered  is 
Gaussian  decay  of  a  Bessel  function: 

p2(r)  =  exp(-R2r2)  J0(Qr),  r_>0  ..  (113) 

Q  equal  to  zero  corresponds  to  Gaussian  spatial  decay  of  the  correlation. 
When  we  substitute  (113)  into  (88),  there  follows  [ll,  6.633  2] 


9This  is  a  form  treated  by  tysanov  (Ref.  17)  and  suggested  by  Schulkin  (Ref.  18,  p.  42). 
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Fig.  4  -  Surface-Height  Spectrum  for  Exponentially  Modulated  Bessel  Function  Spatial  Correlation 
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t 


a\(ij)  =  <t2LxL  UnY 


l  1  /  Q2+<72\f  Qr,\ 

—  «P  ( - - 1 U — r) 

R2  \  4R2  /  \2R2/ 


(114) 


The  surface-height  spectrum  follows  upon  substituting  (114)  into  (92) 
for  the  elliptical  surface  correlation  case.  Since  A^q)  of  ( i 1 4)  is 
positive  for  all  choices  of  Q  and  R,  the  surface-height  spectrum 
4>{f)  of  (92)  is  nonnegative;  therefore,  we  consider  the  spatial  correlation 
form  in  (113)  for  all  Q  and  R. 

The  exponential  and  Bessel  function  in  (114)  can  be  expressed  as 


(115) 


Since  the  fu 'Ction  exp(-t)Io(t)  is  weakly  dependent  on  t,  it  is  seen  that 
Aj  (q  )  possesses  a  peak  approximately  at  q  =  Q  cf  width  proportional 
to  R.  The  integral  of  (92)  for  $(f)  tends  to  smooth  this  peak.  However, 
the  isotropic  surface-height  spectrum  of  (93)  does  not  involve  this 
smoothing  and  is  a  peaked  spectrum;  a  measure  of  the  peakedness  is 
the  "quality"  ratio  of  center  frequency  to  bandwidth,  and  is  related  to 
the  ratio  Q/R  for  this  example. 

The  surface -height  spectrum  for  an  isotropic  surface  is  obtained 
by  substituting  (114)  into  (93a),  and  is  given  by 


where 


x  h  2nf 


■ 

(116) 

£ 

(117) 

(118) 

The  frequency-dependent  term  <|>  is  plotted  versus  the  dimensionless 
parameter  x  in  Fig.  6  for  several  values  of  Q/R.  The  spectrum 
decays  as  exp(-cf4)  for  large  frequencies. 


The  narrow-band  character  of  the  surface-height  spectrum  is 
evident  from  Figs.  4  and  6  when  the  ratio  Q/R  is  large  compared  to 
unity. 
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5  SCATTERING  STRENGTHS  OF  SIDEBAND  COMPONENTS 


For  the  case  of  a  narrow-band  surface-height  spectrum,  the  scat¬ 
tered  acoustic  spectrum  at  a  receiving  point  may  be  obtained  from  (47) 
and  (48).  The  generalform  of  the  received  acoustic  spectrum  was  shown 
in  Fig.  2  to  consist  of  a  series  of  spectral  lobes  or  sidebands,  each 
separated  from  the  incident  acoustic  frequency  fa  by  multiples  of  the 
center  frequency  f  of  the  surface  variation,  plus  an  impulse  at  the 
acoustic  frequency  f a .  This  impulse  is  the  coherent  component. 

The  total  acoustic  power  in  each  spectral  lobe,  or  equivalently  the 
scattering  strength,  may  be  obtained  from  (60)  and  (61)  for  the  spt  cular 
direction,  and  from  (62)  and  (63)  for  the  nonspecular  di-ection.  In 
this  section,  we  will  evaluate  Um(P)  and  Vm(a,|3),  as  given  by  (61) 
and  (63),  respectively.  As  mentioned  in  Subsection  3.4,  if  the  factor 
preceding  U_  ((3)  is  held  constant,  the  geometry,  acoustic  frequency, 
and  correlation  distances  are  considered  fixed.  The  plot  of  Um  (|3), 
therefore,  measures  the  dependence  of  the  mth  scattering  strength 
on  the  surface  rms  wave  height  a  through  the  roughness  parameter 
p.  Vm  (a ,  p )  measures  the  dependence  of  the  mth  scattering  strength 
on  the  correlation  distances  L  and  L  (although  their  ratio  is  fixed) 
through  a,  and  on  the  rms  wave  height  o  through  p.  (Integrals  (61) 
and  (63)  were  numerically  evaluated  by  Simpson's  Rule;  the  error 
analysis  is  contained  in  Appendix  F.) 

In  Figs.  7  through  9,  Um(P)  is  plotted  versus  p  for  m-  0,1, 2,  3,  4 
for  an  exponentially  modulated  Bessel  function  spatial  correlation, 

p2(r)  =  exp(-t)  J0(Qr)  ,  (119) 

for  values  of  Q-  8,4,0.  In  Figs.  10  through  12,  a  similar  set  of  curves 
is  plotted,  the  only  difference  being  that  the  spatial  correlation  is  a 
Gaussianly  modulated  Bessel  function: 

p2( r)  =  e*p(-r2)  Jo(Qr)  ,  (120) 

To  better  explain  these  figures,  let  us  concentrate  temporarily  on 
Fig.  12;  here  Q  =  0,  meaning  the  spatial  correlation  is  Gaussian. 

U,  (P)  corresponds  to  the  scattering  strength  in  the  first  sideband 
(either  above  or  below  the  acoustic  frequency  f4).  Since  a  =  0,  we 
are  considering  only  the  specular  direction,  p  =  0  corresponds  to  a 
flat  surface;  for  this  value  of  P,  there  is  no  scattered  power  in  the 
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Fig.  8  -  Um  for  Exponentially  Modulated  Bessel  Function  Spatial  Correlation,  Q  = 
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Fig.  11  -  U0  for  Gaussianly  Modulated  Bessel  Function  Spatial  Correlation,  Q  =  4 
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sidebands  since  all  the  reflected  power  appears  in  the  coherent  com¬ 
ponent.  As  p  increases,  there  is  less  power  in  the  coherent  component, 
and  more  scattered  power  appears  in  the  sidebands.  The  power  in  the 
first  sideband  reaches  a  maximum  value  when  the  surface  roughness 
parameter  p  is  approximately  unity,  and  decreases  thereafter.  For 
these  larger  values  of  p,  the  surface  is  rougher,  and  more  of  the 
scattered  power  will  appear  in  higher-order  sidebands  in  directions 
other  than  specular.  We  have  plotted  the  zeroth  through  fourth-order 
sideband  scattering  strength  factors  Um(P)  in  Figs.  7  through  12.  All 
the  curves  initially  increase  as  p  increases  from  zero,  reach  a 
maximum,  and  then  decrease.  For  large  p,  there  is  little  difference 
in  power  between  the  various  sidebands. 

For  Q  =  4,  or  8,  the  gross  features  of  the  sideband  powers  are 
similar  to  those  for  Q  =  0,  whether  the  modulation  is  exponential  or 
Gaussian.  However,  whereas  for  Q  =  0,  Uj(P)  is  larger  than  U0(P) 
for  a  large  range  of  p,  the  reverse  is  true  for  Q  =  4  and  8.  For  larger 
Q,  the  zeroth  and  second-order  sidebands  contain  most  of  the  scattered 
power  in  the  specular  direction.  As  Q  increases,  the  power  in  each 
sideband,  for.  a  given  roughness  p,  decreases.  The  reason  for  this 
is  that  the  total  scattered  power  in  the  specular  direction  decreases 
with  increasing  Q,  since  the  correlation  distance  decreases.  (This 
latter  effect  is  due  to  the  factor  J0  (Qr)  in  (119)  and  (120),  which  decays 
with  r  in  addition  to  the  exponential  terms. ) 

For  small  p,  the  power  in  the  first  sideband  is  proportional  to 
p2  ,  and  the  powers  in  the  zeroth  and  second-order  sidebands  behave 
as  p4  (See  (67)).  This  is  not  apparent  in  Fig.  10  for  the  Gaussian 
modulated  case  for  Q  =  8,  because  the  effect  occurs  at  smaller  p 
values  than  plotted. 

The  major  difference  between  the  exponentially  modulated  and 
Gaussianly  modulated  curves  occurs  at  larger  values  of  p.  Here  the 
sideband  powers  are  much  larger  for  the  Gaussian  case  than  for  the 
exponential  case.  This  is  true  for  all  Q  and  all  orders  of  sidebands 
plotted. 

In  Figs.  13  through27,  three-dimensional  plots  of  V  (a,P)  versus 
a  and  p  are  given  for  m  =  0,1, 2, 3, 4  and  for  an  exponentially 
modulated  Bessel  function  spatial  correlation  (See  (119))  for  values  of 
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Fig.  13  -  VQ  for  Exponentially  Modulated  Bessel  Function  Spatial  Correlation,  0  =  8 
(V„  =  .268  •  10" 1  at  a  =  15  1/3,  /3  =  1.6) 
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Fig.  14  -  V  for  Exponentially  Modulated  Besrel  Function  Spatial  Correlation,  0  =  8 
(Vj  =  .763  •  10" 1  at  a  =  8,  /3  =  1.0) 
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Fig.  15  -  V2  for  Exponentially  Modulated  Bessel  Function  Spatial  Correlation,  0  =  8 
(V2  =  .  144  •  10" 1  at  a  =  15  1/3,  /S  =  1.6) 
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Fig.  19  -  V  for  Exponentially  Modulated  Bessel  Function  Spatial  Correlation,  0  =  4 
(V,  =  .395  •  10" 1  at  a  =  4  \fl,  0  =  1.0) 
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Fig.  20  -  V2  for  Exponentially  Modulated  Bessel  Function  Spatial  Correlation,  Q  =  4 

=  .106  •  10*1  at  o-8,  fi  «  1.7) 
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V  for  Exponentially  Modulated  Bessel  Function  bpntial  Correlation,  C  =  •} 


Fig.  24 

(v, 


Vj  for  Exponential  Spatial  Correlation 
.120  •  10* 1  at  a=l  2/i,  fi  =  1.1) 
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Fig.  26  -  V}  for  Exponential  Spatial  Correlation 
(Vj  =  .460  •  10"2  at  a  =  14  1/3,  0  =  3.2) 
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Fig.  27  -  V4  for  Exponential  Spatial  Correlation 


Q  =  8,4,0.  In  Figs.  28  through  42,  a  similar  set  of  curves  is  plotted; 
the  only  difference  is  that  the  spatial  correlation  is  a  Gaussianly 
modulated  Bessel  function  of  the  form  of  (120).  These  functions  were 
computed  for  a  =  0(1/3)20  and  p  =  0(.  1)6.  The  maximum* values 
attained  by  Vm  (a,|3)  at  this  grid  structure  are  presented  in  Table  1 
and  are  noted  on  each  plot.  (These  values  are  not  necessarily  the 
maximum  values  of  Vm  (a,  (3),  but  simply  the  maximum  values  at  the 
sample  points  investigated.)  Since  the  plots  are  isometric,  values  of 
Vm  (a,P)  at  other  values  of  a,  P  can  be  obtained  by  measuring  the 
vertical  distance  above  the  a,  (3  plane  and  scaling  according  to  the 
values  given. 

To  explain  the  detailed  behavior  of  these  plots,  it  is  helpful  to 
consider  (3  as  a  measure  of  the  surface  vertical  roughness,  and  a 
as  a  measure  of  the  surface  horizontal  roughness.  As  (3  increases, 
the  surface  becomes  rougher  in  the  vertical  direction;  also  as  a 
decreases,  i.  e. ,  the  correlation  distances  decrease,  the  surface 
becomes  rougher  in  the  horizontal  direction.  For  a  given  a,  the 
behavior  of  (a,  (3)  with  (3  is  similar  to  that  observed  in  Figs.  7 
through  12. 

For  (3  =  0,  Va  (a,  (3)  is  zero  for  all  a  since  the  surface  is  flat 
and  there  is  no  scattered  power.  As  (3  increases,  the  scattered  power 
increases,  reaching  a  maximum.  As  (3  increases  still  further,  the 
curve  decreases,  indicating  that  more  of  the  scattered  power  appears 
in  higher  order  sidebands.  Note  from  the  plots  that  for  a  given  sur¬ 
face  spatial  correlation  and  a  given  Q,  the  location  (a,  (3),  at  which 
Vm(a,P)  reaches  a  maximum  value,  increases  as  m  increases;  i.  e. . 
the  peak  location  recedes  from  the  origin. 

For  a  given  (3  (fixed  roughness)  and  a  large  a  (large  correlation 
distances),  the  surface  is  very  planar,  resulting  in  most  of  the  re¬ 
flected  power  appearing  in  the  specular  direction.  As  a  decreases 
from  large  values,  the  surface  becomes  less  planar,  more  power  is 
scattered  in  the  nonspecular  direction,  and  Va  increases.  However, 
as  a  decreases  still  further,  the  surface  becomes  rougher  in  the 
horizontal  direction,  and  more  of  the  scattered  power  appears  in  the 
higher  order  sidebands,  so  that  each  particular  Vm(a,P)  decreases 
with  a  in  this  range.  The  ratio  p/a  is  a  measure  of  the  average 
slope  of  the  surface;  as  this  ratio  increases,  more  of  the  scattered 
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Table  1 


APPROXIMATE  LOCATIONS  OF  MAXIMUM 
VALUES  OF  V,  (q,(3) 


Fig. 

a 

No. 

13 

15  I  '3 

14 

8 

15 

15  1/3 

16 

20 

17 

20 

18 

7  2'5 

19 

4  1/3 

20 

8 

21 

17 

22 

20 

23 

4 

24 

1  2  3 

25 

6  2  3 

26 

14  1  3 

27 

20 

28 

15 

29 

8  1/3 

30 

15  1/3 

31 

20 

32 

20 

33 

7  2/3 

34 

4  2/3 

35 

8 

36 

13  i'J 

37 

18 

38 

3 

39 

2 

40 

4 

41 

6  1/3 

42 

8 

(q,P) 


.2  68 
.763 
.  144 
.678 
.392 
.  196 
.395 
.  106 
.601 
.425 
.  106 
.  120 
.643 
.460 
.352 
.325 
.705 
.  177 
.890 
.570 
.249 
.398 
.  141 
.903 
.690 
.210 
.237 
.  125 
.886 
.683 


10* 1 
10* 1 
10- 1 
lO*2 

io-2 
10- 1 
10-* 
10- 1 
10~2 
io-2 
io-1 
io-1 
io-2 
10-2 
10-2 
10-» 
io-1 
io-1 

IO*2 

IO*2 

IO*1 

io-1 

IO-1 

IO*2 

io'2 

IO"1 

10** 

IO-1 

10-2 

10-2 


Fig.  28  -  Vq  for  Gaussiaulv  Modulated  3essel  Function  Spatial  Correlation,  0  =  8 
(V„  -  .325  •  IQ- 1  at  a  »  15,  /3  -  1.6) 
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Fig.  29  -  V,  for  Gaussianly  Modulated  Bessel  Function  Spatial  Correlation,  Q  =  8 
(Vt  =* .705  •  HT*  at  a  *8  i/i,  ft  m  1.0) 
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unction  Spatial  < 


30  -  Vj  fot  Gaimianly  Modulated  Bessel  i- 


F;g.  31  -  Vj  for  Gaussianly  Modulated  Bessel  Function  Spatial  Correlation,  0  =  8 
(V3  =  .890  ■  10“2  at  a  =  2G,/3~  2.7) 
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i<  M'  JfcJu^MkWiVl  A  A'l  Mi  >'> 


Fig.  32  -  V4  for  Gautsianly  Modulated  Be<vl  Function  Spatial  Correlation,  0  =  8 
(V4  *  .570  •  10* 2  at  a  =  20,  /3  =  3-l) 
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Fig,  34  -  V,  for  Gauaianly  Modulated  Bessel  Function  Spatial  Correlation,  Q 
(y  j  =  .398  •  10" 1  at  a  *  4  2/3,  /3  =  l.l) 


Fig.  36  -  V3  for  Gaussianiy  Modulated  Bessel  Function  Spatial  Correlation,  0  =  4, 
(Vj  =  .90?  .  10' 2  at  a  =  13  i/j,  fi  =  3- 1) 
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Fig.  42  -  for  Gaussian  Spatial  Correlation 
(V4  =  .683  •  10”2  at  a  e  8,  /3  -  4.1) 
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power  appears  in  the  higher  order  sidebands.  This  behavior  can  occur 
through  either  increased  £  (roughness)  or  decreased  a  (correlation 
distances). 

As  a  approaches  zero,  for  the  nonspecular  direction,  the  corre¬ 
lation  distances  approach  zero.  For  this  case,  the  surface  becomes 
very  rough  and,  in  the  limit  of  zero-correlation  distances,  has  the 
properties  of  white  noise.  Although  mathematically,  for  a  =  0,  the 
amount  of  scattered  energy  in  the  sidebands  is  zero,  the  initial  as¬ 
sumptions  (See  Section  l)are  violated  for  small  correlation  distances; 
therefore,  this  portion  of  the  plots  does  not  have  physical  significance. 

Curves  of  Vm  (a,|3)  for  p2  (r)  =  exp(-r)  cos(Qr),  Q  =  1,  were  also 
obtained,  but  have  not  been  included  because  they  are  very  similar  to 
the  corresponding  curves  for  p2  (r)  =  exp(-r)  J  (Qr);  in  particular,  the 
omitted  curves  are  more  peaked  than  exp(-r),  but  not  as  peaked  as 
exp(-r)  Jo  (4r).  Similarly,  curves  of  Vm  (a,£)  for  p2  (r)  =  exp(-r2)* 
cos(Qr),  Q  =  1.  8482777,10  were  also  obtained  and  exhibited  behavior 
intermediate  to  that  of  exp(-r2)  and  exp(-r2 )  Jo  (4r). 

Let  us  now  compare  the  exponentially  modulated  Bessel  function 
plots  with  those  of  the  Gaussian  modulated  Bessel  function  plots.  For 
large  £  ,  the  Gaussian  case  has  more  power  in  a  given  sideband  than 
that  of  the  corresponding  exponential  case.  This  behavior  is  similar 
to  that  of  the  curves  in  Figs.  7  through  12.  For  Q  =  0,  the  corre¬ 

sponding  exponential  and  Gaussian  plots  are  somewhat  dissimilar. 
However,  for  Q  =  4,  the  corresponding  curves  of  the  exponential  case 
are  more  nearly  similar  to  those  of  the  Gaussian  case,  whereas  for 
Q  =  8,  thetwosets  of  curves  are  almost  identical.  Thus,  for  small  Q, 
the  envelope  of  the  spatial  correlation  of  the  heights  is  the  controlling 
factor,  butas  Q  increases,  the  number  of  oscillations  per  given  distance 
of  Jo  (Qr)  increases,  and  Jn  (Qr)  becomes  the  controlling  factor  in 
spatial  correlation. 


6.  DISCUSSION 

In  summary,  an  analysis  was  made  of  a  time-varying  random 
surface  statistically  stationary  in  both  space  and  time.  The  reflected 

,0Tkese  choices  of  Q  are  die  largest  possible  consistent  with  a  nonnegative  surface- height  spectrum; 
see  Subsection  4.3. 
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acoustic  spectrum  was  obtained  in  terms  o£  the  first-  and  second -order 
characteristic  functions  of  the  surface -height  variation.  For  the  special 
case  of  a  Gaussian  distribution  of  surface  heights,  a  narrow-band 
surface  spectrum,  and  a  surface  spatial-temporal  correlation  function 
that  is  stationary  in  the  wide  sense,  the  mathematical  equation  for  the 
received  acoustic  spectrum  showed  that  the  scatter  spectrum  is  com¬ 
posed  of  spectral  lobes,  or  sidebands,  centered  at  frequencies  equal 
to  the  acoustic  frequency  plus  multiples  of  the  surface -center  frequency 
(See  (47)  and  (48)).  For  the  special  case  of  a  surface  spatial  correla¬ 
tion  that  has  elliptical  contours,  the  evaluation  of  the  complete  spectrum 
can  be  expressed  as  a  Fourier  transform  of  a  single  integral.  The  com¬ 
plete  spectrum  was  not  numerically  evaluated  in  this  report,  but  should 
be  obtained  in  future  work.  In  addition,  the  restriction  to  a  symmetric 
directional  wave  spectrum,  (82),  should  be  eliminated.  For  example, 
received  spectra  for  hemispherical  directional  wave  spectra  should  be 
computed;  this  will  yield  unequal  scattered  sideband  powers. 

The  power  in  each  sideband  was  evaluatedfor  a  variety  of  conditions 
and  spatial  correlations.  It  was  found  that  the  zeroth  and  second-order 
sideband  powers  had  a  very  similar  behavior  with  surface  roughness, 
whereas  the  other  sidebands  had  a  somewhat  different  dependence.  It 
has  been  assumed  that  there  is  not  an  appreciable  overlap  of  the  side¬ 
bands;  for  large  surface  roughness,  there  would  be  an  appreciable 
overlap  of  sidelobes,  and  the  complete  spectrum  would  have  to  be 
obtained  for  this  case. 

It  is  worthwhile  to  discuss  the  difference  between  the  surface 
treated  in  this  report  and  the  fixed -amplitude  sinusoidal  surface  con¬ 
sidered  by  Roderick  and  Cron  [7j.  We  have  considered  a  narrow-band 
spectrum  for  the  height  variation  at  apointon  the  surface,  and  assumed 
the  joint  probability  density  of  the  surface  heights  to  be  Gaussian.  As 
the  bandwidth  of  the  surface  variation  decreases  and  approaches  zero, 
the  properties  of  this  surface  process  do  not  approach  the  fixed-am¬ 
plitude  sinusoidal  surface  case.  For  example.,  the  distribution  of 
heights  remains  Gaussian  and  is,  therefore,  different  from  the  prob¬ 
ability  density  associated  with  a  sinusoid.  The  scattered  powers  in 
the  sidebands,  as  given  by  the  narrow- band  Gaussian  theory,  are 
different  from  those  of  the  single-frequency  sinusoidal  theory. 
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In  mathematical  terms,  the  first-order-characteristic  function  of  a 
fixed-amplitude  sinusoid  is  a  Bessel  function.  If  we  average  this  Bessel 
function  with  the  Rayleigh  distribution  of  sine  wave  amplitudes,  we 
obtain  the  Gaussian  first-order-characteristic  function.  Thus,  our  re¬ 
sults  can  be  viewed  as  a  particular  average  over  the  sinusoidal  surface 
case. 

However,  a  narrow-band  time  function  with  center  frequency  f 
and  bandwidth  W  resembles  a  sine  wave  of  frequency  ft  over  a  short 
period  of  time,  but  will  differ  from  a  sine  wave  over  a  longs’*  period 
of  time.  The  difference  will  show  up  on  the  order  of  W_1  sec.  Thus, 
if  we  obtain  a  short  sample  of  the  reflected  signal  from  a  narrow-band 
surface,  and  determine  the  received  acoustic  spectrum  from  it,  the 
single-frequency  theory  (with  the  current  surface -height  amplitude) 
should  accurately  predict  that  spectrum.  However,  for  a  long-time 
average  over  the  statistical  properties  of  the  surface,  the  narrow -band 
Gaussian  theory,  not  the  fixed -amplitude  single -frequency  theory,  must 
be  used  to  predict  the  expected  behavior  of  the  spectrum. 

This  report  has  concentrated  on  the  spectral  quantities  of  the  re¬ 
flected  signal.  However,  the  probability  distributions  of  the  envelope 
of  the  individual  sidebands  can  also  be  obtained  in  at  least  one  important 
case.  If  the  area  of  insonification  is  much  larger  than  the  correlation 
area  of  the  surface  (product  of  correlation  distances),  we  can  appeal 
to  the  Central  Limit  Theorem  and  state  that  the  envelope  and  phase  dis¬ 
tributions  of  the  mth  component  of  the  received  scatter  pressure  field 
are  Rayleigh  and  uniform,  respectively,  for  m  /  0.  For  m  =  0,  the 
scatter  component  combined  with  the  coherent  component  results  in  a 
Beckmann  distribution  [2,  Chapter  7j.  (Of  course,  as  the  coherent 
component  tends  to  zero,  the  Beckmann  distribution  tends  to  the  Rayleigh 
distribution.)  It  should  be  noted  that  the  assumption  of  Gaussian  surface 
statistics,  i.  e. ,  £{x,  y,  t)  Gaussian,  does  not  necessarily  imply  that 
the  received  pressure  is  Gaussian,  because  the  received  pressure  is 
obtained  via  a  nonlinear  transformation  on  C  (See  (3));  rather,  it  is 
the  assumption  of  a  large  number  of  independently  reflecting  surface 
portions  that  yields  the  Gaussian  behavior. 
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Appendix  A 
GEOMETRY  FACTOR 


The  factor  B  in  (3)  is  given  different  values  by  different  authors. 
For  example,  the  values  for  B  given  by  ParkinsAI  and  Beckmann  and 
Spizzichino,A2  using  the  notation  in  this  report,  are 


1  +  uft 

B  = - SL 


“a  _  1+aQaA  +  bQbA  +  CQCA 


CQ  +  CA 


whereas  EckartAJ  gives  the  following  value  for  B: 


(A-l) 


B  = 


CQ  +  CA 


(A-2) 


Beckmann  and  Spizzichino*4  attribute  the  discrepancy  to  the  replace¬ 
ment  of  the  normal  to  the  surface  by  a  vector  in  the  z-direction,  an 
approximation  that  is  valid  only  for  surfaces  with  very  gentle  slopes. 


For  the  specular  direction, 

aA  =  “aQ*  bA=“bQ’  c A  =CQ 


(A -3) 


both  forms  for  B  above  become  cQ.  Thus,  there  is  no  difference  in 
the  specular  direction.  However,  there  can  be  significant  differences 
for  other  cases.  For  example,  the  second  form  for  B  depends  only  on 
the  direction  cosines  with  respect  to  the  z-axis,  but  the  first  form 
depends  on  all  direction  cosines.  Thus,  at  normal  incidence. 


aQ  =  bQ=0’cQ  =  1 


(A -4) 


a1B.  E.  Parkins,  "Scattering  from  the  Tune- Varying  Surface  of  the  Ocean, "  Journal  of  the 
Acoustical  Society  of  America,  vol.  42,  no.  6,  December  1967,  Eq.  (10),  p.  1263. 

Beckmann  and  A.  Spizzichino,  The  Scattering  of  Electromagnetic  Waves  from  Rough 
Surfaces,  The  Macmillan  Co.,  New  York,  1963,  Eqs.  (11)  and  (12),  p.  27,  and  Eqs.  (29)  and  (32), 
pp.  22-23. 

A5C.  Eckart,  "The  Scattering  of  Sound  from  the  Sea  Surface, "  Journal  of  the  Acoustical 
Society  of  America,  vol.  26,  no.  3,  May  1953,  Eq.  (6),  p.  567. 

A4 

Beckmann  and  Spizzichino,  Appendix  A,  op.cit. 
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and  the  two  forms  for  B  become 


l  and  '/$(l  +  cA)  ,  (A-5) 

respectively.  For  an  observation  point  near  grazing,  these  terms  differ 
by  a  factor  of  2  in  pressure. 

Rather  than  attempt  to  resolve  the  differences  in  the  scale  factor, 
we  have  expressed  all  our  quantities  in  terms  of  the  general  symbol  B 
and  leave  it  to  the  reader  to  make  his  choice  as  to  the  correct  factor. 
For  a  fixed  geometry,  so  far  as  the  relative  spectral  content  of  the 
received  acoustic  waveform  is  concerned,  the  exact  geometric  factor 
is  not  important  under  the  assumptions  of  the  present  theory. 


Appendix  5 

CORRELATION  PROPERTIES  OF  SINGLE-SIDED  PROCESS 


Suppose  x(t)  is  a  wide-sense  stationary  single-sided81  (complex) 
process  with  no  dc  component.  Let 


jt{t)  =  xf(t)  +  ix.(t) 


(B-l) 


where 

*tw  =  a  f*w +«•(*)]  . 

Jtj(t)  =  ^{x(t)-x*(t)l 


(B-2) 


Denoting  the  Hilbert  transform  by  U I  I ,  we  have 

HU.WUKtHfcWl +  #}«•«» 

=  %[-i  x  (t)  +  i  x*  (t)l  =  x.  (t) 


(B-3) 


where  the  single-sided  character  of  x(t)  and  its  lack  of  dc  has  been 
utilized.  The  Hilbert  transform  is  represented  by  a  linear  network  with 
transfer  function  -isgn(f).  Therefore,  the  transformation  between 
xt(t)  and  x;  (t)  can  be  represented  as  in  Fig.  B-l. 


«r  ft) 


-itgn(f) 


*i  ft) 


Fig.  B-l  -  Hilbert  Transform 

Thus,  we  have  derived  the  fact  that  the  imaginary  part  of  a  wide-sense 
stationary  single -sided  process  is  the  Hilbert  transform  of  the  real 
part.®2  When  the  properties  of  linear  networks  are  used,  it  follows 
from  Fig.  B-l  that 


81  Single-tided  meant  that  the  power  density  spectrum  it  confined  to  positive  frequencies. 

®2  This  is  the  converse  of  the  usual  situation  where  a  waveform  plus  i  times  its  Hilben  transform 
is  shown  to  be  tingle-tided. 


81 


x((t)  x.(t-r)  =  Xf(t)  xf(t-r)  , 

xf(t)  s  (t-r)  =  -x.(t)  xf  (t-r) 

Combining  these  properties,  we  find  that 

x(t)  x(t-r)  =  0  for  all  r  . 

Of  course,  x(t)  x*  (t-r)  /  0, 
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Appendix  C 

SCATTERING  COEFFICIENT  AND  SCATTERING  STRENGTH 


The  received  (complex)  pressure  p(t)  at  the  observation  point  was 
given  in  (3).  The  scatter  component  ps(t)  was  defined  in  (11),  and  the 
mean-square  value  of  p^  (t)  was  evaluated  in  (18)  and  (20).  In  this 
appendix,  we  will  relate  the  scattering  coefficient  and  scattering  strength 
to  these  quantities. 


The  average  scatter  intensity  (power/unit  area)  at  the  receiver  for 
a  real  sinusoidal  signal  transmittal  is 


I 


S 


Ip,U>I2/2 

pv 


(C-l) 


where  p  is  the  density  of  the  fluid  and  v  is  the  propagation  velocity.01 
If  the  receiver  at  A  subtends  a  solid  angle  Q  ,  as  seen  from  origin  O 
in  Fig.  1,  the  average  received  scatter  power  is 

IsR20  ,  (C-2) 

and  the  average  received  scatter  power  per  steradian  is 

,*R»=^7lp*(t)l2Ro  ’  (C'3) 

spy 

Since  the  incident  pressure  on  the  surface  at  x,y  is  p.  (x,  y) 
(See  (3)),  the  incident  intensity  at  x,  y  is 

.  (C-4) 

2  pv 

Therefore,  the  total  incident  power  on  the  surface  is 

cQ //dxdy— i-|p.(x,y)  |2  =— 2-J  (0,0)  ,  (C-5) 

w  2pm  1  2  pm 

using  (16). 


The  scattering  coefficient  or  scattering  cross  section  as  is  defined 
by  Eckart02  as 


Cl  The  factor  of  i/2  follows  in  a  manner  similar  to  that  in  (5)  and  (6). 

02 C.  Eckart,  "The  Scattering  of  Sound  from  the  Sea  Surface,"  journal  of  the  Acoustical 
Society  of  Am .•rica,  vol.  25,  no.  3,  May  1953,  Eq.  (18),  p.  568,  and  Eq.  (10),  p.  567. 


83 


*s=' 


2  d  2 
o 


i  P.  <«>  1 2  R, 
J  (0,0) 

which  can  be  manipulated  into  the  forms 


(C-6) 


a  = 


|p,(t) 


I  2  o  2 


I  R 

%  O 


/jfd*dy  Ipj(x.y)  |2  ff  dx  dy  _i_  (x,  y)  | 2 

Ipv 

average  received  scatter  power  per  steradtao 


(C-7) 


incident  power  on  surface  /  cQ  ’ 

where  we  have  used  (16),  (C-3),  and  (C-5). 

For  purposes  of  this  report,  the  "scattering  strength"  S  is  defined 

I 


as 


C3 


(C-8) 


where  I  t  is  the  average  scatter  intensity  at  the  receiver  due  to  unit 
scattering  area,  referred  to  unit  distance,  and  I.oc  is  the  average  inci¬ 
dent  intensity  on  the  surface.  If  we  let  Ae(|  be  the  effective  area  of 
insonification  in  the  x,  y  plane  (e.  g. , 


l*K 


//dxdy 


P=  (0,0) 


)■ 


then 


using  (C-l).  Also, 


2  q2 

“o 


i  bJ  Ip.WI2  R, 

— s t  R* 

K  0 
ycfl 


2pr  Aeff 


(C-9) 


(C-10) 


incident  power  on  die  surface 


i»e 


cQ  J(0,0) 


(C-ll) 


2 pV  Aelf 

using  (C-5).  Taking  the  ratio  of  (C-10)  and  (C-ll),  we  note  that  the 
scattering  strength  becomes 


S  = 


IPS(*>I2  Ro  _ff» 
c Q  J(0,0)  “«Q 


(C-12) 


CJ  R.  J.  Urick,  Principles  of  Underwater  Sound  for  Engineers,  McGraw-Hill  Book  Co.,  New  York, 
1967,  p.  108,  Section  8.2.  Actually,  Urick  has  a  decibel  ratio  rather  than  a  linear  power  ratio., 
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using  (C-6).  Thus,  the  scattering  strength  equals  the  scattering  coef¬ 
ficient  divided  by  the  direction  cosine  of  the  incident  pressure  with  re¬ 
spect  to  the  normal  to  the  surface.  This  latter  quantity,  os  /cQ,  which 
is  exactly  what  has  been  recommended  by  some  authors04  as  prefer¬ 
able  to  a t  alone,  has  been  adopted  in  this  report. 


Notice  that  the  exact  definition  of  Ae(f  didnotenter  into  S  because 
it  cancelled  in  the  ratio  (C- 8);  therefore,  (C-9)  is  merely  one  possible 
definition  of  this  quantity.  An  alternate  definition  of  S,  which  avoids 
this  quantity  altogether,  is  available  upon  substitution  of  (C-7)  into(C-lZ): 


S 


average  received  scatter  power  pet  stcradian 
incident  power  on  surface 


(C-13) 


«C.  W.  Horton,  Sr.,  and  T.  G.  Muir,  "Theoretical  Studies  of  the  Scattering  Acoustic  Waves 
from  a  Rough  Surface, "  journal  of  the  Acoustical  Society  of  America,  vol.  41,  no.  3,  March  1967, 
Section  U.C.,  p.  630. 
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Appendix  D 

COMPLEX  ENVELOFE  REPRESENTATION  OF 
NARROW- BAND,  SURFACE-HEIGHT-CORRELATION  FUNCTION 


From  (40),  the  surface -height  correlation  is 

a2  p( u,  v,  r)  =  £(x,y,  t)  £(x- u,  y-v,  t-r) 


(D-l) 


These  are  all  real  quantities.  When  the  homogeneity  and  stationarity 
of  the  surface  is  utilized,  there  follows 

p(-u, -v,  — r)  =  p(u,  v,  t)  .  (D-2) 

We  define  a  normalized  cross -spectrum  g  as  the  Fourier  transform 

°*  P"  g(u,v,f)  =  /drexp(-i2nYr)  p(u,v,  r)  ^ 

It  is  a  function  of  position  differences  u,  v,  as  well  as  frequency  f, 
and  is  neither  real  nor  even  in  f;  however,  it  is  narrow-band.  Let  us 
now  define  the  single-sided  spectrum  as 

/  A  <2g(u, v,  f),  f>0l 


t  a  <2g(u,v,f),  f>0l 

o, ( <o)  • 


and  the  center  frequency  of  surface  variations  at  a  point  as 


/dff  g+(0,0,f)  £dffg(0,0,f) 

£  _  .  ,  .  /  _  '  X  —  bo 


(D-4) 


(D-5) 


f  dfg(0,0,f) 


The  expression  g(0,  0,  f)  is  proportional  to  the  surface -height  spectrum 
at  a  point;  it  is  real  and  even  because  p(0,  0,r)  is  real  and  even.  We 
now  shift  the  single-sided  spectrum  g+  down  by  fs  to  get  a  low -pass 
spectrum  centered  about  zero  frequency: 

j:<u,v,f)  =  g+(u,v,f+fj)  .  (D-6) 

Then,  the  Fourier  transform  of  g  is  defined  as  the  complex  envelope  (&: 

(o(u,  v,  r)  =  / df  exp (i 2»r  f  r) _g(u,  v,  r) 

=  p+(u,v,  r)  e*p(-i2irftr)  ,  ^  ^ 

where  we  have  used  (D-6)  and  defined  p+  as  the  Fourier  transform  of 
g  +  .  Now,  using  f*  Fourier  transform  of  the  unit  step,01 


f  df  exp  (i2irft)  =  lA  S(t)  +  - 


(D-8) 


0lA.  Papouiis,  The  Fourier  Integral  and  Its  Applications,  McGraw-Hill  Book  Co.,  New  York, 
1962,  Eq.  (3- 13). 
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(D-9) 


p+  can  be  evaluated  <”om  (D-4)  as 

p+(u,v,  r)  =  p(u,v,  r)  +  ipH(u,  v,  r)  , 

where  pH  is  the  Hilbert  transform  of  p: 

1  \  1  £  p(u,v,w) 

pH  (u,  v,  r)  =—  fdw -  . 

H  nj  r-  w 

Combining  (D-9)  and  (D-7),  we  obtain 

p(u,  v,  r)  =  Re{p+(u,  v,  r)} 

=  Rela(u,  v,  r)  exp(i2fffB  r)l  , 
which  is  the  property  desired  in  (44). 


(D-10) 

(D- 11) 


An  additional  property  of  complex  envelope  a  that  will  be  necessary 
is  now  derived  from  (D-10): 

pH(-u,-v,-r)  =-pH(u,v,  r)  .  (D-12) 

Substituting  into  (D-9),  we  note  that 

p+(-u,-v,  -r)  =  p*  (u,  v,  r)  ,  (D- 13) 

and  then  from  (D-7),  we  get 

g(-u,-v,-r)  «g?(u,v,  r)  .  (D-14) 

For  the  special  case  where  p(u,v,r)  isevenin  r  for  all  u,  v,  i.  e. , 

p(«,T,~r)=p(o,v,  r)  ,  (D-I5) 

ft(u,  v,  0)  is  purely  real  and  is  given  by 

£(o,t,0)  =p(u,v,0)  .  (D-16) 

This  may  be  seen  as  follows:  From  (D-10),  using  the  evenness  of  p, 

pH(o,v,0)  *0.  (D-17) 

Then,  using  (D-9)  and  (D-7),  (D-16)  follows.  , 
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Appendix  E 

NUMERICAL  EVALUATION  OF  EQ.  (112) 


In  (112),  the  surface -height  spectrum  was  given  in  terms  of  an 
integral: 

$U)  »  |x|3  /  dt  t  exp(-t2)  cos(qt)  J0(*2t)  ,  (E-l) 

0 

where  q  a  Q/R.  To  evaluate  this  integral,  we  expand  cos(qt)  in  apower 
series  in  t  and  obtain 

3U)  =  M3BIo(-l)n  /  *«2o+Ie^>(-t*)JoU2t)  .  (E-2) 

The  integral  in  (E-2),  from  Gradshteyn  and  Ryzhik,61  is 

Hn!«p(-*V4)Ln(K4/4)  ,  <E-3> 

where  Ln  is  a  Laguerre  polynomial.  Ther  .ore, 

$(*)»H|*|3exp(-*4/4)  2  (— q2)”  ~~L  (*4/i)  .  (E-4) 

>°o  (2nj! 

This  series,  which  can  be  easily  evaluated  by  using  recurrence  relations 
on  the  Laguerre  polynomials,  is  plotted  in  Fig.  5. 


For  x  near  zero,  $(x)  becomes  negative  if  q  is  made  too  large. 
From  (E-l),  it  may  be  seen  that  the  largest  value  of  q  allowed  is  when 

f  it  t  expt-t2)  cos  (qr)  *  0  .  (E-  5) 

That  is,  from  Gradshteyn  and  Ryzhik,EJ 

M  jFjtiuK;  -Jfq2)  c0  .  (E-6) 

The  solution  for  the  smallest  value  of  q  satisfying  (E-6)  is 

q  1.8482777  .  (E-7) 

Therefore,  (112)  has  been  plotted  only  for  Q/R  less  than  this  value. 


E1I.  S.  Gradshteyn  and  I.  W.  Ryzhik,  Table  of  Integral!,  Series,  and  Products,  Academic  Pte», 
New  York,  1965,  Eq.  6.631  10. 

E2  Ibid.,  Eq.  3.952  8. 


89/90 

REVERSE  BLANK 


Appendix  F 

ERROR  ANALYSIS  OF  EQS.  (61)  AND  (63) 


The  infinite  integral  of  (61)  must  be  approximated  by  a  finite  integral, 
up  to  value  L.,  say.  In  order  to  determine  how  large  L  must  be  for 
negligible  error,  we  note  that  for  small  x 


I  to  -S  ~ 

tn  on 


(x/2)2  ,  m  =0 


(x/2)“ 

—  ,  m.>  1 


On!)2 


(F-l) 


Therefore,  for  large  values  of  r. 


I„02Patt>-*o.~ 


B4/32p2(r))2  ,  m=0 


(F-2) 


The  slowest  decaying  order  is  m  =  1;  therefore,  the  largest  error  in 
terminating  the  integral  (61)  at  L  is  approximately 

E=(2*r)“1e«p  (— /S2)  /  drr%/32p.(r)  .  (F-3) 

L 

When  spatial  correlation  p2  is  an  exponential  modulation  of  either  a 
cosine  or  a  Bessel  function, 

p2(r)  =  exp(-r)  cos(Qr)  or  e*p(-r)  Jo(Qr)  ,  (F-4) 

(F-3)  is  upper-bounded  by  M 

E  £.(2ff)-1  expi-p2)  MfS2  /  dr  r  exp(-r) 

L  (F-5) 

*=  (4ff)“ 1  P2  exp  (-/32)  (L  + 1)  exp  (-L)  . 

The  largest  value  of  p2  exp(-p2)  occurs  at  P  =  1,  yielding 

E  <.(4ffe)_l  (L  +  l)  exp(-L)  . 

For  L  =  25,  E  <.I0"U.  The  limit  25  was  used  in  evaluating  (61)  for 
the  exponentially  modulated  spatial  correlation;  all  the  values  of  Um  in 
Figs.  7,  8,  and  9  are  much  larger  than  10' ll. 


For  a  Gaussianly  modulated  spatial  correlation, 

p2(r)  =  exp  (-r2)  cos  (Qr)  or  exp  (-r2)  Jo  (Qr)  ; 


(F-7) 
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the  bound  corresponding  to  (F-6)  is 

E  <_(8ffe>-1  e*p(-L2)  .  (F-8) 

For  L  =  5,  E<2  x  10-13.  The  limit  of  5  was  used  in  evaluating  (61) 
for  the  Gaussianly  modulated  spatial  correlation;  all  the  values  of  Um 
in  Figs.  10,  11,  and  12  are  much  larger  than  2  x  10”13. 

For  integral  (63),  we  can  upper-bound  the  integrand  by  replacing 
the  Bessel  function  by  unity.  The  slowest  decaying  order  is  again  m  =  1. 
An  analysis  similar  to  that  in  (F-l)  through (F-3)  yields,  for  exponential 
modulation,  the  upper  bound  on  the  error 

E  <(4ffe)_1  a2(L+l)  e*p(-L)  ,  (**-9) 

whereas,  for  Gaussian  modulation,  the  upper  bound  on  the  error  is  given  by 

E  staffer1  a2  eipt-L2)  .  (F-10) 

For  L  =  25,  (F-9)  yields  10"ua2;  for  L  =  5,  (F-10)  yields  2  x  10~«a2. 

The  interval  (0,  L),  used  in  evaluating  (61)  and  (63),  was  repeatedly 
cut  in  half  and  evaluated  by  Simpson's  rule  until  insignificant  change 
occurred  in  the  sum. 
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The  power  density  spectrum  of  the  received  acoustic  waveform  after  reflection  off 
the  time-varying  random  sea  surface  is  evaluated  at  an  arbitrary  observation  point  in 
die  farfield.  For  a  monochromatic  transmitted  signal  and  a  narrow -band  Gaussian 
surface -height  variation,  the  received  acoustic  spectrum  is  shown  to  consist  of  an 
impulse  at  the  transmitted  acoustic  frequency  plus  sidelobes  centered  at  frequencies 
separated  from  the  transmitted  frequency  by  multiples  of  the  surface  center  frequency. 
The  powers  in  the  coherent  component  and  scattered  sideband  components  of  the 
received  pressure  waveform  are  evaluated  in  terms  of  the  surface  roughness  and 
spatial- temporal  correlation  function  of  the  surface.  For  the  special  case  of  elliptical 
contours  of  iso-correlation  at  zero  time  delay,  the  sideband  powers  and  scattering 
strengths  are  evaluated  in  terms  of  two  fundamental  parameters  that  include  the 
geometry  of  the  experiment,  the  incident  acoustic  frequency,  the  root  mean  square 
(rms)  surface  height,  and  the  surface  correlation  distances. 

The  rms  bandwidths  of  the  sideband  scatter  components  are  evaluated  for  small 
surface  roughness  and  shown  to  be  approximately  proportional  to  the  square  root  of 
the  sideband  number.  Numerous  examples  of  sideband  scattering  strengths  for  a 
variety  of  spatial  correlation  functions,  including  exponential  and  Gaussian  decay  as 
special  cases,  are  given. 
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